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Abstract 



We study the potential of a Neutrino Factory in constraining the parameter space 
of a scheme with one sterile neutrino separated from three active ones by an 0(1) 
eV^ mass-squared difference. We present approximated analytic expressions for the 
oscillation probabilities, showing that the greatest sensitivity to sterile neutrinos at 
a Neutrino Factory can be achieved using the — > and the z/^ — > v-r oscillations. 
We have studied two setups: a Neutrino Factory with 50 GeV (20 GeV) stored 
muons, with two detectors of the Hybrid-MIND type (a magnetized ECC next to 
a magnetized iron calorimeter), located at L = 3000,7500 km (L = 4000,7500 
km) from the source. Four channels have been used: Ue i^^iji^r', t^^i i^fi^i^r- 
The relevant backgrounds, efficiencies and systematic errors have been taken into 
account, and we have discussed dependence of the sensitivities on the systematic 
errors. We have found that the 50 GeV (20 GeV) setup can constrain sin^ 29^^^'^'^^ < 
7x10-^(2x10"'^); (934 < 12°(14°); and 624 < 7.5°(8°). Our results hold for any value 
of Am^BL ^0.1 eV^ Eventually we have shown that, if a positive signal is found, the 
proposed setup is able to measure simultaneously ^34 and S3 with a precision of few 
degrees and few tens of degrees, respectively, solving the so-called "intrinsic" and 
"sign degeneracies". Combination of disappearance and of the — > Uj- channel, 
that will be called " the discovery channer , at the two baselines is able to measure 
at 99% CL a new CP-violating phase 63 for sin^ 2^34 > 0.06. 



1 Introduction 

From the results of solar pH2|3|^f5|6|71l8] . atmospheric [QllTO] . reactor pT|12|13l[n] and 
accelerator [T5]ll6lll7p!8] neutrino experiments, we now know that neutrinos have masses 
and mixings. In the framework of three flavor oscillations, neutrino oscillations are de- 
scribed by three mixing angles, 6^12, 6^13, 623, and one CP phase 6, as well as two inde- 
pendent mass-squared differences, Arrigj^ and Amgj^. In the standard parametrization 
jl9j for the three flavor mixing matrix t/pMNS pilll21|22|[23j . when 613 is small, {Am^i, 
612) and (Amgj^, 6*23) correspond to the mass-squared difference and the mixing angle 
of the solar and atmospheric oscillations, respectively. From the solar neutrino exper- 
iments we have Am^^i ~ 7.7 x 10"^ eV^ sin^ 6112 ~ 0.30 - 0.31 [21125] . and from the 
atmospheric neutrino experiments |Am^^j^| ~ 2.4 x 10"^ eV^, sin^ 2^23 — 0.47 — 0.50 
[21125] . As for ^13, the reactor data [TT] [12fl3j and three-family global analysis of the 
experimental data give an upper bounJJ~sii?6'i3 < 0.04, while we have no information 
on 5 at present (see, however, Ref. [28]). 

In order to determine precisely the remaining two parameters ^13 and 5, long baseline 

^ In Refs. |26|25|27] . a global analysis of the neutrino oscillation data has been shown, in 
which a non-vanishing value for ^13 is found. This result, however, is compatible with ^13 = 
at less than 2a, and it has not been confirmed by the other groups performing global fits [24j . 
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experiments with intense neutrino beams have been proposed |29ll30f31ll32ll33f34ll35j . As 
in the case of the B factories these precision measurements will allow us to look 

for deviation from the standard three flavor oscillations scenario. Possible scenarios for 
such deviations include non-standard interactions which affects the neutrino produc- 
tions and detections [38], those which modify the neutrino propagations [39|^ . light 
sterile neutrinos [22], unitarity violation due to the effect of heavy fermions [11112], 
etc. These scenarios (except the non-standard interactions which change the neutrino 
propagation, only), in general break unitarity of the PMNS matrix. As in the B physics, 
test of unitarity is one of the important problems which should be investigated in the 
future long baseline experiments (see Ref. |13] for a review). Among the proposed long 
baseline experiments with high intense neutrino beams, a Neutrino Factory [34], which 
uses a muon storage ring to produce neutrino beams from muon decays, is expected to 
have excellent sensitivity to ^13 and S. 

One of the advantages of a Neutrino Factory is that the flux is flavor-rich, well under 
control and with no contamination. This last point is of particular relevance for new 
physics searches in neutrino oscillations, since the Ur oscillations provide one of 

the most promising signal of non-standard physics in oscillations (as it will be shown in 
Sec. 12.41 for sterile neutrinos; see also Refs. [i^45f46|H7] ) and z/^ detection is important 
to check unitarity of the PMNS matrix (although, clearly, unitarity violations of the 
PMNS matrix are best studied in weak decay processes, @Ij). A Neutrino Factory with 
multi-GeV muons is a powerful facility to look for r's signals, if detectors dedicated 
to r-detection are provided. Notice that oscillations into z/^ cannot be measured by 
so-called (3 beams [35J and that high energy conventional superbeams are affected by 
Pt- contamination of the flux (through 5-mesons decay). 

Four-neutrino mass schemes have attracted much attention since the announcement by 
the LSND group on evidence for neutrino oscillations with a mass squared dif- 

ference Am^ ~ 0(1) eV^ |l8]l49f50j . Because the mass squared difference suggested by 
the LSND result is much larger than those for the solar and atmospheric neutrino oscil- 
lations, in order to explain all these data in terms of neutrino oscillations, it is necessary 
to introduce at least a fourth light neutrino state. From the LEP data [T9|5T] . which 
indicate that the number of weakly interacting light neutrinos is three, the fourth state 
has to be a sterile neutrino. For this reason, the LSND signal could be considered as an 
evidence for the existence of a sterile neutrino. Recently the MiniBooNE experiment 
|52j gave a negative result for neutrino oscillations with the mass squared difference 
Am^ ~ 0(1) eV^ which was suggested by the LSND data, and it has become difficult 
for four-neutrino models to explain the LSND data. The so-called (3+2)-scheme with 
two sterile neutrinos has also been proposed [5^ to account for LSND, but also in 
this case, tension with the disappearance experiments remains, as long as we take into 
account the LSND data. Adding a third sterile neutrino does not seem to help [51] , 
and in general global analyses seem to indicate that sterile neutrinos alone are not 
enough to account for all the data in terms of neutrino oscillations. Models with sterile 
neutrinos and exotic physics have been therefore proposed [5^56ll57|58|59j . 
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While the efforts to account for all the data including the LSND in terms of neutrino 
oscillations have been unsuccessful, sterile neutrino scenarios which satisfy all the ex- 
perimental constraints except LSND are still possible. Even if the inclusion of light 
sterile neutrinos is not needed to explain the present experimental data, it is certainly 
worth investigating scenarios where sterile neutrinos do appear and constrain their pa- 
rameter space. Light singlet fermions are indeed present in the low-energy spectrum of 
many theories and models including them represent, for example, phenomenologically 
natural frameworks to break three-flavor unitarity. 

In Ref. ^Q] the (3-1-1 )-scheme without imposing the LSND constraint was studied in the 
context of the CNGS experiments ^j, finding that if the OPERA detector is exposed 
to the nominal CNGS beam intensity, a null result can improve a bit the present bound 
on ^13, but not those on the active-sterile mixing angles, 614,624 and 6^4. 

In this paper, we have extended the analysis in Ref. [60] to the case of a Neutrino 
Factory experiment. We have first of all extended the analytic computation of the 
oscillation probabilities for the (3 -f l)-model at long baseline experiments in matter 
using the formalism by Kimura-Takamura-Yokomakura (KTY) [B^IIUH] . Approximated 
formulae in powers of ^13, of the deviations from maximality of 623 {S623) and of the 
active-sterile mixing angles, 6i4, have been obtained. On the basis of this analysis, we 
have found that the greatest sensitivity to the active-sterile mixing angles is achieved 
using the z/^ — > z/^ and —>■ v^- channels (as it was noticed, for example, in Refs. [MfUH] 
and refs. therein). To take full advantage of these signals, detectors capable of both z/^ 
and Vr identification are needed. In our numerical analysis we have, thus, assumed a 
detector of the Hybrid-MIND type [66]: a 50 kton magnetized iron calorimeter next to 
a 4 kton Emulsion Cloud Chamber with magnetized iron plates. This detector has a 
greater efficiency to than the standard OPERA-type ECC, with lead plates 

acting as target. 

We have then extensively analyzed the physics reach of a 50 GeV Neutrino Factory 
which has 2 x lO^'' useful muon decays per year aimed at two detectors of the Hybrid- 
MIND type located at L = 3000 km and L = 7500 km from the source, with both 
polarities running for 4 years each. As a consistency check, we have also studied the 
case of a 20 GeV Neutrino Factory which has 5 x lO^'' useful muon decays per year 
aimed at the same two detectors of the Hybrid-MIND type located at L = 4000 km and 
L = 7500 km, again with both polarities running for 4 years each. The latter option 
was the scenario which was suggested in the International Scoping Study for a future 
Neutrino Factory and Super-Beam facility [13] . 

Four signals have been considered: the "standard" Neutrino Factory channels, the 
golden channel z/g [67] and the silver channel Ve [68j; the disappear- 

ance channel; and the novel signal u^-, that will be named in this paper the 

"discovery channeT . We think that this name is appropriate due to its high-sensitivity 
to new physics signals in neutrino oscillations, as found from the theoretical analysis of 
oscillation probabilities. At the same time, however, we are forced to remind that the 
effective potential of this channel is tightly linked to the performances of the proposed 
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detectors able to identify r's, that are still in a very preliminary design phase and that 
should be studied further. 

Using the first two channels at the 50 GeV Neutrino Factory, we can extend the three- 
family ^13-sensitivity plots to the four-family (6'i3^^™'', 6'i4)-plane, putting a stringent 
upper bound on 6[^^^^ mixing angle, sin^ 2^^3^^™'' < O(10~*^) at 90% CL. Using the 
combination of the z/^ disappearance channel and of the " discovery channel " z/^ —>■ Vr 
at the 50 GeV Neutrino Factory, we are able to constrain the active-sterile mixing 
angles ^34 and ^24, \0m\ < 12° and |^24| < 7.5° at 90% CL. We have found that the 
combination of the shortest baseline data with the longest baseline ones significantly 
improves the sensitivity in both planes, {0^ii'^'^\Oii) and {92A,9^i). 

We have, then, compared the results for the 50 GeV Neutrino Factory with those that 
can be obtained at the 20 GeV ISS-inspired one (with 2 x 10^° and 5 x 10^° useful 
muon decays per year per baseline, respectively). We have found that the former setup 
has a greater potential than the latter for sterile neutrino searches, in particular for 
the simultaneous measurement of 62^ and 6*34. We interpret this result as consequence 
of the larger r statistics that can be collected using higher energy muons. 

We have also studied the region of the four-family parameter space for which a four- 
neutrino signal cannot be confused with the three-family model, determining the "dis- 
covery potential" of the Neutrino Factory in the [Oi'^ \,9ii)- and (6*24, 6'34)-planes. 

We have, then, performed a preliminary study of the potential of the 50 GeV Neutrino 
Factory to measure four-family parameters, showing 99% CL contours in the (^34, ^3)- 
and (^13^'^™'*, 52)-planes. We have found that the shortest baseline data are affected by 
"intrinsic" [69] and "sign" degeneracies [70]. On the other hand, the longest baseline 
data have a rather good precision in ^34 and are not affected by those degeneracies. 
The combination of the two baselines, thus, is able to measure simultaneously ^34 and 
(^3 with a precision of a few degrees and a few tens of degrees, respectively. In this 
context, we have also studied the " (53-discovery potential". We have found that the 
disappearance channel can measure a CP-violating ^3 for sin^ 2^34 > 0.4 with a 50% 
CP-coverage; the combination of the disappearance and of the z/^ "discovery" 
channels, however, is able to measure a CP-violating ^3 for for sin^ 26'34 > 0.06 with a 
80% CP-coverage. The combination of the two channels is extremely effective in solving 
correlations and degeneracies, thus improving the 53-discovery potential by an order of 
magnitude in sin^ 26'34. We have also shown how the measurement of the three-family 
like CP-violating phase is modified by the presence of non-vanishing active-sterile 
mixing angles. 

Eventually, we have studied the impact of statistic and systematic errors on the per- 
formance of the detector proposed to identify r's (the magnetized ECC component of 
the Hybrid-MIND detector, MECC). We have found that, to improve significantly the 
potential of the MECC, systematic errors should be kept at the level of a few percents. 
On the other hand, increasing the mass of the MECC from 4 kton to 8 kton is found 
to be of marginal impact. It is clear from this analysis that the r-detector section of 
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(a) (b) 

Fig. 1. The two classes of four-neutrino mass spectra, (a): (2+2) and (b): (3+1). 

the Hybrid-MIND should be studied further to understand its ultimate potential. To 
this purpose, we will take great advantage of the understanding of the ECC technology 
and systematics after the first years of OPERA data taking. 

An analysis such as this is not completely new: first studies of sterile neutrinos at a 
Neutrino Factory were presented at the first NuFact workshop in Lyon in 1999 [7T|72] 
in the framework of the so-called (2+2)-schemes and subsequently extended to the case 
of (3+l)-schemes in Ref. The possibility to use the Neutrino Factory detectors, 
optimized to look for three-family oscillations, to disentangle three- from four-neutrino 
signals was considered in Ref. [73)165] . Recently, in Ref. ^74j a four-family neutrino 
analysis in the spirit of Ref. [60] has been performed. The main differences between 
this paper and Ref. [71] are (i) that careful numerical analyses are carried out here 
by taking into account backgrounds, efficiencies and systematic errors specific to the 
considered signals and setup, and (ii) that the four channels i/^ — >• z/,-, z/^ — z/^, — >• z/^, 
z/g at the Neutrino Factory are considered and their contributions are clarified in 
the present paper. 

The paper is organized as follows. In Sec. [2] the main features of four-neutrino schemes 
and the present bounds on the mixing angles in these scenarios are briefiy summa- 
rized. Furthermore we compute approximated oscillation probabilities in matter in the 
atmospheric regime using the KTY formalism [5^|U3] (details of our computations are 
given in App. A). In Sec. [3] we remind the details of the considered Neutrino Factory 
setup. In Sec. Hlwe present our results using various channels at the Neutrino Factory. 
Finally, in Sec. [5] we draw our conclusions. 



2 Four neutrino schemes 

Four-neutrino schemes consist of one extra sterile state in addition to the three weakly 
interacting ones. Depending on whether one or two mass eigenstate(s) are separated 



6 



from the others by the largest mass- squared gapifj, the schemes are called (3+1)- and 
(2-|-2)-schemes, as is shown in Fig. [TJ In the (3+1) schemes, there is a group of three 
close- by neutrino masses that is separated from the fourth one by the larger gap. In 
(2+2) schemes, there are two pairs of close masses separated by the large gap. These 
two classes lead to very different phenomenological consequences. 

2.1 (2+2)-schemes 

A characteristic feature of (2+2) schemes is that the extra sterile state cannot si- 
multaneously decouple from both solar and atmospheric oscillations. The fraction of 
sterile neutrino contributions to solar and atmospheric oscillations is given by rjs = 
\Usi\^ + \ Us2\^ and l — rjs = |f^s3p + |f^s4p, respectively, where the mass squared differ- 
ences Am2i and |Am43| are assumed to be those of the solar and atmospheric oscilla- 
tions. The experimental results show that mixing among active neutrinos give domi- 
nant contributions to both the solar and atmospheric oscillations (see, e.g., Ref. [7S]). 
In particular, in Fig. 19 of Ref. [76] we can see that at the 99% level r]s < 0.25 and 
1 — rjs < 0.25, which contradicts the unitarity condition I]j=i Wsjl"^ = 1- In fact the 
(2+2)-schemes are excluded at 5. la CL [76]. This conclusion is independent of whether 
we take the LSND data into consideration or not and we will not consider (2+2)- 
schemes in the rest of this paper. 

2.2 (3+l)-schemes with the LSND constraint 

On the other hand, (3+l)-schemes are not affected by the tension between the solar and 
atmospheric constraints on sterile neutrino oscillations, because as long as the mixing 
of sterile neutrino is small, then phenomenology of solar and atmospheric oscillations 
is approximately the same as that of the three flavor framework. The (3+1) schemes 
start having a problem only when one tries to account for LSND and all other negative 
results of the short baseline experiments. To explain the LSND data while satisfying the 
constraints from other disappearance experiments, the oscillation probabilities of the 
appearance and disappearance channels have to satisfy the following relation [77|75] : 

sin^ 26'LSND(Am^) < ^ sin^ 26'Bugey(Am^) ■ sin^ 26'cDHsw(Am^) (1) 

where ^LSND(Am^), ^cDHsw(Am^), 6'Bugey(Am^) are the value of the effective two-flavor 
mixing angle as a function of the mass squared difference Am^ in the allowed region for 
LSND (z/^ — i> z/e), the CDHSW experiment [79] (i/^ Ufj), and the Bugey experiment 

^ The only assumption for the largest mass-squared difference is that oscillations caused 
by this mass-squared difference are averaged. So the results hold for any value of Am^gL ^ 
0.1 eV^. Interesting models with "sterile" neutrino with masses m ~ 0(1) KeV can be found, 
for example, in Ref. |75j . 
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[80] (z/g — ^ Ue), respectively. The reason that the (3+l)-scheme to explain LSND is 
disfavored is basically because eq. ([T]) is not satisfied for any value of Am^. A (3+2)- 
scheme with two sterile neutrino has also been proposed [S^ to account for LSND, and 
it may be possible to reconcile the LSND and MiniBooNE data by introducing a CP 
phase [8T|M] . Also in this case, however, tension with CDHSW [79] and Bugey [HO] 
remains, as in the case of the (3+l)-scheme. 

2.3 (3+l)-schemes without the LSND constraint 

If we give up our effort to account for the LSND data, on the other hand, we no 
longer have the constraint ([I]). In this case we have only the upper bound on the extra 
mixing angles and this scenario satisfies all the experimental constraints (except that 
of LSND). Throughout this work, therefore, we will consider a (3+l)-scheme without 
taking the LSND data into account while satisfying all the negative constraints, as it 
was done in Ref. [HQ] . 

It has been discussed that the mixing angles of four neutrino schemes may be con- 
strained by big-bang nucleosynthesis (see Refs. [77|82] and references therein), and if 
such arguments are applied, then the mixing angles of sterile neutrinos would have 
to be very small. However, it is known that in some model [HS] neutrino oscillations 
themselves create large lepton asymmetries which prevent sterile neutrinos from being 
in thermal equilibrium, so it is not so clear whether the arguments in [771IH2] hold. At 
present, therefore, it is fare to say that there is not yet general consensus on this issue 
(see Ref. [HI] and references therein). In this paper we will not impose cosmological 
constraints on our scheme. 

The mixing matrix U can be conveniently parametrized in terms of six independent 
rotation angles 6ij and three (if neutrinos are Dirac fermions) or six (if neutrinos 
are Majorana fermions) phases 5j. In oscillation experiments, only the so-called "Dirac 
phases" can be measured, since the "Majorana phases" appear only as an overall phase 
of the oscillation amplitude and disappear in the oscillation probability. The Majorana 
or Dirac nature of neutrinos can thus be tested only in AL = 2 transitions such as 
neutrino- less double /5-decay or lepton number violating decays [I9]. In the following 
analysis, with no loss in generality, we will restrict ourselves to the case of 4 Dirac-type 
neutrinos only. 

A generic rotation in a four-dimensional space can be obtained by performing six 
different rotations along the Euler axes. Since the ordering of the rotation matrices 
Rij (where ij refers to the plane in which the rotation takes place) is arbitrary, plenty 
of different parametrizations of the mixing matrix U are allowed. In this paper we 
are interested in the so-called "atmospheric regime", with oscillations driven by the 
atmospheric mass difference, Aatm = ~ 0(1). We will then make use of 

the following parametrization, adopted in Ref. [M] : 

U = RuiOu, 0) i?24(^24, 0) i?23(^23, ^s) i?14(^14, 0) RrsiOis, 62) i?12(^12, ^l) • (2) 



8 



In eq. ([2]), Rij{9ij, Si) are the complex rotation matrices in the ij-plane defined as: 

cos 6ij P = Q = h j 
1 p = q^ij 

[RijiOij, = \ sin 9i, e-''^ P = i]q = ] (3) 

— sin 6ij e^^' p = j]q = i 
otherwise. 

It is convenient to put phases in R12 (so that it automatically drops in the limit Asoi = 
Am'^^iL/2E — > 0) and R13 (so that it reduces to the "standard" three- family Dirac 
phase when sterile neutrinos are decoupled). The third phase can be put anywhere; 
we will place it in R23. Note that in the one-mass dominance regime [S3] (i.e. for 
Aatm, Asoi 0) all the phases automatically disappear from the oscillation probabilities 
for some choices of the four-family PMNS matrix parametrization. 

The mixing matrix elements in the parametrization ([2]) are given in the Appendix El 

This parametrization has been used in Ref. |60j to put bounds on the active-sterile 
mixing angles 6^ using existing data (including MiniBooNE but neglecting LSND). 
These bounds can be summarized as follows: 

(1) Bounds from disappearance reactor experiments 

Reactor experiments such as Bugey and Chooz can put stringent bounds on 
6'i3 and 614 in this parametrization: 613 < 13° and 614 < 10° at 99% CL, for any 
value of Am^gL > 0.1 eV^, with some correlation between the two (in particular 
the four-family Chooz bound on 6*13 is slightly modulated by 6*14). 

(2) Bounds from disappearance experiments 

A "negative" result in a z/^ disappearance experiment at "atmospheric" L/E 
(such as K2K or MINOS), in which z/^ oscillations can be very well fitted in terms 
of three- family oscillations, puts a stringent bound on the mixing angle 624. The 
bound from such experiments found in Ref. [60] is: 624 < 14° at 99% CL, for 

Am^^L > 0.1 eV^ 

(3) Bounds on 634 

Neither z/g nor z/^ disappearance probabilities in vacuum depend strongly on 634 
(as it can be seen in Ref. [nO])- An upper bound on 634, however, can be drawn 
as the result of indirect searches for Vg conversion in atmospheric neutrino 

experiments, that take advantage of the different interaction with matter of active 
and sterile neutrinos. Present bounds on 634 arise, thus, from a measurement of 
spectral distortion. On the other hand, bounds on ^^13, 614 and ^24 are mainly drawn 
by a flux normalization measurement. As a consequence, the bound on 634 that we 
can draw by non-observation of z/^ — > Vg oscillation in atmospheric experiments 
is less stringent than those we have shown before. For this reason, 634 can be 
somewhat larger than 613,614 and 624- 634 < 32° at 99% CL. 
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These bounds are depicted in Fig. 2 of Ref. |60], where 90%, 95%, 99% and 3(t CL 



contours in the (6*13 — O14)- and (6^24 — ^34)-planes are shown for Asoi and Aml^^ 
2.4 X 10-3 eV^ and Oo^ = 45°. 



2.4 Oscillation probabilities at the Neutrino Factory in the (3+l)-scheme 

To understand the details of the different channels with the greatest sensitivity to the 
four-family neutrino schemes, it is useful to obtain simple analytical expressions for the 
different channels in matter. Hereafter, we will assume a constant Earth density along 
the neutrino path, computed using the PREM |86]. Notice that, in the framework of 
Neutrino Factory experiments simulations, the impact of non-constant matter density 
has been thoroughly studied, showing that for the baselines and the muon energy 
considered in this paper the details of the density profile crossed by neutrinos does not 
modify the resultsjiJ 

To derive the oscillation probabilities, we adopt the KTY formalism [62)163] (the details 
are given in the Appendix [B]) . Furthermore, to get simplified forms of the formula, it 
is convenient to obtain the probabilities expanding with respect to the following small 
parameters: 



^'34 ~ v^'ia ~ V^i4 ~ v6'24 ~ \J 5(^2z < 4 X 10 ^ 



with 5^23 = ^23 — 7r/4. 

The relevant oscillation probabilities in matter, expanded to third order in e, are given 

by 

P,e~l + 0(e^), (4) 

Pe^^Per ~ ^e. ~ O (e^) , (5) 



Pt,^, = l- sin^ ~ 2 (AiL) S24 S34 cos sin A31L + O (e"^) , (6) 

P^r = (1 - ^34) sin^ + {s24 S34 sin S3 + 2 S24 S34 cos ^3} sin A31L 

+0 {e') , (7) 

= sLsin^ — ^ - ■524S34sin(53sin A31L + O (e"^) , (8) 



where A31 = Am3^/2£', and we take the convention that the central value of [Amg^l 
is Am^jj^, which is determined by the two flavor analysis of the atmospheric neutrino 



3 This is totally different, for example, in the case of /3-Beams with Li/B decaying ions. In 
those setups, a resonance is crossed for 0(5 GeV) neutrinos for L = 0(10000) km baselines 
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data. The matter density parameter An is An = y^Gpn-n/^. Notice that at O(e^) the 
expansion parameter 5^23 is not present in the oscillation probabilities (it only arises 
at the next order in e). 

From eqs. (Hj)-®, it can be easily verified that unitarity of the PMNS matrix is satisfied 
to this order in e. As it can be seen, moreover, the u^. decouples within this approxi- 
mation. We can thus conclude that the "classic" Neutrino Factory channels, such as 
the "golden channel" —>■ and the "silver channel" z/g — > Vri are of limited inter- 
est to study sterile neutrinos, as we will see later in the numerical analysis!^ Leading 
sensitivity to 6*34 is provided by the first term in P^^, eq. ([7]), that is proportional to 
(1 — s1^). Sensitivity to ^24 is best achieved using the z/^ — > z/^ oscillations, though 
the relevant term appears at O(e^), as it will be shown in Sec. 14.21 In eqs. dH])-© we 
can also see that combined sensitivity to 624, and 6^4^ is achievable through the O(e^) 
matter-dependent S24S34COS53 term in Pf^^.P^j^r- The bounds on these two angles are 
those that can be improved the most by the Neutrino Factory experiments. We can 
thus safely say that the disappearance channel and the — »• Vr appearance channel 
are the most relevant signals to look for sterile neutrinos. This will be confirmed by 
the numerical analysis later. 

Notice that ^34 is considerably less constrained than the rest of the new angles (which 
are as constrained as ^13). Therefore, the z/^ — v-r channel would have at the theoretical 
level the strongest sensitivity to sterile neutrino parameter space: P^r is the only one 
(with the exception of P^g, that cannot be directly measured) that is of O(e^) in the 
expansion parameters. Moreover, z/^ — »• z/^ can also look for CP- violating signals in four- 
family scenarios through the O(e^) matter-independent S24S34sin53 term, something 
out of the reach of the z/^ disappearance channel (that can only measure ^3 through 
CP-conserving signals). 

As it is well known, r-detection experiments are extremely difficult at the experimental 
level. If the experimental problems could be overcome, would be the important 

channel to study sterile neutrinos as well as other kinds of new physics, such as the 
non-standard interactions [^^11^ . or unitarity violations [lUllTF] . In particular, if CP 
violation occurs due to these new physics, then the — > channel is quite powerful 
in measuring the new CP-violating phases. For this reason, we will name it throughout 
this paper as " the discovery channeT . 

Notice that the — > v-r channel has not been studied in detail yet in the framework 
of neutrino factories studies (with the possible exceptions of Refs. [88|89] ). 

3 The Neutrino Factory and the Hybrid-MIND detector 

In a Neutrino Factory [Mf90] muons are first produced with a multi-MW proton source, 
accelerated up to energies of several GeV's, and finally injected into a storage ring 

^ This has been known since long. See for example Refs. |64|65|60] . 



11 




with long straight sections aiming to one ore more detectors. The muon decays jj^ — 
e"*" Ue Ufj_ and fi^ e~ Ue v^i provide a very well known two-flavor neutrino flux [91] 
with energies in the range Ey e [0,£'/i]. Neutrino Factory designs have been proposed 
in Europe the US (HUESinSIHZ] , and Japan [98j. The dedicated International 

Scoping Study for a future Neutrino Factory and Super-Beam facility [43j showed that, 
provided suflicient resources, an accelerator complex capable of providing about 10^^ 
muon decays of a given polarity per year can be built. 

The Neutrino Factory setup that we propose for sterile neutrinos searches and that 
we examine in detail in the rest of the paper is deflned as follows: muons of both 
polarities are accelerated up to = 50 GeV and injected into one storage ring with a 
geometry that allows to aim at two far detectors, the flrst located at 3000 km and the 
second at 7500 km from the source. An alternative option, considered in the flnal ISS 
Accelerator Report [99], is to inject the muon beam into different storage rings, each 
of them aimed to a single far detector. The number of useful muon decays per year 
aimed at each detector has been flxed to 2 x 10^°. This number is rather conservative 
since, in the flnal ISS Physics Report [13], with a similar storage ring(s) geometry, 
5 X 10^° useful muon decays per year aimed at each detector are considered (i.e., 10^^ 
total useful muon decays per year). Four years of data taking for each muon polarity 
are envisaged. 

The Neutrino Factory fluxes (for fi~ accumulated in the storage ring) at L = 3000 km 
as a function of the neutrino energy for E^ = 50 GeV are shown in Fig. [2]|^left). 

Two detectors of different technologies have been considered to detect the i/^ and Ur 
signals. The flrst one is a magnetized iron calorimeter [102] . that was proposed with 
a slightly different design in Ref. [103] to measure the "golden" z/g — >■ i'^ wrong-sign 
muons signal. The second detector is a Magnetized Emulsion Cloud Chamber (MECC) 
[66], an evolution of the ECC modeled after OPERA that was flrst considered for 
Neutrino Factory studies in Refs. [68|104] to look for the "silver channel" z/g —>■ Vt (a 
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channel that can be looked at only in a Neutrino Factory setup). 

The v^N and Vt-N cross-sections on iron as a function of the neutrino energy for both 
neutrinos and anti- neutrinos are shown in Fig. [2I^right). 

Notice that the adopted setup, similar to those proposed in the first Neutrino Fac- 
tory studies (see, for example, Refs. [571IU^ ). slightly differs from the setup that was 
suggested in the final ISS Physics Report. The latter consists in stored muons with 
energies in the range e [20, 30] GeV, aiming to two detectors located at L = 4000 
km and L = 7500 km from the source. The longest baseline corresponds to the so-called 
"magic baseline" |105j . where three-family CP- violating effects vanish. We have chosen 
the same baseline for the far detector in our setup, although in the (3-1-1 )-neutrino 
model not all of the CP-violating phases dependence decouple at this distance (as it 
can be seen from eqs. ( 16)171) and in Sec. H]). The shortest baseline was optimized in the 
ISS Physics Report to look for CP- violating three-family signals, finding that a detector 
with a baseline of L = 4000 km performed slightly better than at L = 3000 km. This 
optimization, however, is no longer valid when looking for the (3+l)-neutrino model 
signals. We decided, therefore, to adopt the L = 3000 km baseline used in previous 
studies, for which possible sites have already been exploredlZI Eventually, we use a 
stored muon energy that is larger than the optimal value adopted in the ISS Physics 
Report, again chosen to maximize the sensitivity to three- family observables such as 
^13, the sign of the atmospheric mass difference and the three-family CP- violating phase 
S. It is indeed well known from previous studies (see, for example, Ref. |il06j for an op- 
timization of the muon energy to look for NSI signals at the Neutrino Factory) that to 
look for new physics the higher the muon energy the better. An evident motivation for 
this is that the z/^ — » z/^ channel is very important to look for new physics in neutrino 
oscillations, and high neutrino energies are required to circumvent the extremely low 
i>T-N cross-section in the tens of GeV energy range (see Fig. [21 right). 

To show that the 50 GeV Neutrino Factory setup proposed above is more suited to look 
for sterile neutrino signals, we will compare our results with an ISS-inspired Neutrino 
Factory, defined as follows: muons of both polarities are accelerated up to = 20 
GeV and injected into storage ring(s) with a geometry that allows to aim at two far 
detectors, the first located at 4000 km and the second at 7500 km from the source. 
The number of useful muon decays per year aimed at each detector has been fixed in 
this case to 5 x 10^°, following the final ISS Physics Report [43j. Four years of data 
taking for each muon polarity are envisaged. When studying the performance of the 
ISS-inspired 20 GeV setup, we consider the same detectors as in the case of the 50 GeV 
Neutrino Factory. 



^ For example, the island of La Palma in the Canary Islands, the island of Longyearbyen in 
Norway or the Oulu mine in Pyhosalmi in Finland are possible options for the location of a 
detector at approximately 3000 km from the source for a CERN-based Neutrino Factory. 
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3. 1 The MIND detector: z/g and z/^ 

The baseline detector, a 50 kton magnetized iron calorimeter of the MINOS type, 
was originally optimized to reduce the background to the wrong-sing muon signal for 
T^e i^fi oscillations (represented dominantly by right-sign muons with a wrong charge 
assignment and charmed meson decays) to the 10~^ level. To achieve this extremely 
ambitious signal-to-background ratio, tight kinematical cuts were applied. Such cuts, 
although strongly reducing the background, have the disadvantage of an important 
suppression of the signal below 10 GeV, an energy region that, on the other hand, has 
been shown to be extremely important. A non- negligible signal below and above the first 
oscillation peak (that, for L ~ 3000 km, lies precisely in this energy range) is crucial 
to solve many of the parametric degeneracies [69]I70|107II108] that bother the three- 
family (6'i3, 5)-measurement. To reduce such problems, a modification of the original 
detector proposal called MIND (Magnetized Iron Neutrino Detector) was presented in 
Ref. [102j . 

Studies of the four-family Uf, oscillation (still the best channel to measure 9i^) 

are not different from those performed in the framework of the three-family model. 
In particular, no new sources of background are expected. Quite the contrary, those 
backgrounds induced by wrongly reconstructed i/^ are expected to decrease for large 
values of 6*34, due to the increased oscillation into sterile neutrinos, see eq. ((H]). When 
looking at oscillations, we will therefore take advantage of the wrong-sign muon 

identification efficiency presented in the ISS Detector Report [66j: e^^ = 0.7 above 10 
GeV, with the efficiency increasing linearly from eefi = 0.1 at 1 GcV. 

The MIND detector can also be used to look for z/^ — > z/^ disappearance (one of the 
channels of interest to study sterile neutrino models), providing a very good measure- 
ment of the atmospheric parameters 623 and Aml^^ and giving some handle to solve the 
"octant degeneracy" (see, e.g., Ref. [SS|109] ) . For the right-sign muon sample there 
is no need to accurately tell the charge of the muon, since the background induced 
by misidentified wrong-signs muons is negligible with respect to the signal, [llOj . We 
can safely use for this signal the muon identification efficiency of the MINOS experi- 
ment [111] : e^^ = 0.9 above 1 GeV. Notice that at the MIND detector it is not possible 
to single out r's decaying to muons. We cannot thus use MIND to study the leading 
i^fi oscillation and the "silver" z/g — ^ z/^ channel. 

We have considered as background for the disappearance channel 10~^ of all neutral 
current events, all wrong-sign muon events and the right-sign muons coming from 
oscillation with r decaying into muons. The inclusion of this background has 
no effect on our results for this channel, that are remarkably systematic-dominated. 

Different treatments of the energy response of the detector can be found. For example, 
in Ref. [6?] a constant energy resolution l\Ey = 0.2E^ was considered, by grouping 
events in five bins of energy width AE^ = 10 GeV. On the other hand, in Ref. |112] a 
finer binning was adopted, with a more refined treatment of the energy resolution: 43 
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bins of variable AE^, were considered in the energy range G [IGeV, E^], folding the 
event distribution with a GauBian resolution kernel of variable width, aE = 0.15 x E,^. 
In this paper we have followed the first approach, grouping events into 10 constant 
energy resolution bins with AEi, = O.li?^, leaving possible improvements of the detector 
simulation following Ref. |112j (if needed) for future publications. 

Throughout our numerical simulations we have assumed 2% and 5% for the bin-to- 
bin uncorrelated systematic errors on the golden channel and on the i/^ disappearance 
signal. We have also assumed 1% and 5% for normalization and energy spectrum dis- 
torsion as the correlated systematic errors for all the channels!^ The dependence of 
the sensitivities on the systematic errors will be discussed in section 14.41 

Notice that a different proposal for a magnetized iron calorimeter that can be used for 
a Neutrino Factory experiment has been advanced in Ref. |113] . 

3.2 The MECC detector: u^. — >• Ur and — > Ur 

Two technologies were considered in the literature to study neutrino oscillations into 
r's: Liquid Argon detectors [88f89] and Emulsion Cloud Chamber techniques. In both 
cases, the z/g, z/^ signal can be tagged looking for right-sign muons in coincidence 
with a r decay vertex, to distinguish them from disappearance muons. Therefore, a 
detector with muon charge identification and vertex reconstruction is needed. 

A dedicated analysis to use an ECC modeled after OPERA at the Neutrino Factory 
to look for the "silver channel" z/g [[68] has been published in Ref. |104j . In that 

reference, a 5 kton ECC was considered, with a detailed study of the main sources of 
background. 

The result of that analysis was that r's can be identified with a very low background, 
at the price of a very low efficiency, 0(5%). One of the main motivation of such low 
efficiency was that only the r — >■ /i decay channel was used, i.e. only 17% of the total 
amount of produced r's. In the OPERA detector used at the CNGS and in the analysis 
of Ref. |104] . the r charge identification is achieved using two large spectrometers 
located at the end of two thick active sections, with bricks made of repeated layers of 
lead (acting as the target for v^N interactions) and emulsions. In the Magnetized ECC 
(MECC) proposal [66J the lead plates are replaced by iron plates, again interleaved 
with emulsions layers. The ECC can be, thus, directly magnetized through the iron 
plates. Emulsion spectrometers (currently in their test phase, |114j ) are located at 
the end of the ECC section. Eventually, the MECC is placed in front of the MIND 
detector, to form the so-called Hybrid-MIND setup. The efficiency of this detector to 

^ In principle, some of the correlated systematic errors could be common among the different 
channels (such as the detector volume uncertainty) or different baselines (such as the cross 
section uncertainty), but for simplicity we assume here that all the correlated systematic 
errors are independent among the different channels or different baselines. The results in the 
present paper may be, therefore, somewhat conservative. 
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r's is much higher than in the case of the standard ECC, since r decay into electrons 
and into hadrons can be used in addition to r — /i decays. The expected efficiency is, 
thus, approximately five times larger than in the case of the ECC (see, for example, 
Refs. [6611115] and refs. therein). 




The MECC bricks are bigger than the corresponding ECC ones, due to the replacement 
of the lead target with iron. The huge volume to be magnetized put, thus, a tight limit 
on the maximum foreseeable detector mass. We will therefore consider throughout the 
paper a 4 kton MECC |116j located in front of the 50 kton MIND, and use this detector 
to study both the silver channel and the novel "discovery channel" — >■ v^-. 

For the silver channel signal, we will use an energy dependent efficiency eer taken from 
Ref. |104] . multiplying it by a factor five to take into account that all r's decay channel 
can be used at the MECC. A detailed study of the efficiency to the z/^ — > v-r channel 
at the MECC, on the other hand, is lacking. We will therefore assume a constant 
efficiency e^,- = 0.65 above 5 GeV, increasing by a factor five the efficiency considered 
in Ref. [10] • This assumption must be checked in further studies of the MECC- type 
detectors exposed to a Neutrino Factory beam. 

The backgrounds for the silver and the discovery channels should be also correspond- 
ingly increased at the MECC with respect to the ECC ones. At the ECC, the expected 
signal-to-background ratio (after some kinematical cuts) for (using the t —>■ fi 

decay channel, only) is 50:1 or larger [UU], the dominant source of background for 
the process z/^ — > i/^ — > r~ — >• fi~ being represented by non-oscillated muons that 
produce charmed mesons eventually decaying into /i~ either through NC or CC in 
which the muon is not observed. No detailed study of the expected background for the 
i/g Ur or the — > Ur signals at the MECC exposed to a Neutrino Factory beam 
has been performed yet, though. We have thus decided to make the assumption that, 
using MECC, all r decay channels should be affected by similar backgrounds. We have 
therefore consistently multiplied the backgrounds for z/g — > and z/^ zA_computed 
in Refs. |104|60] by a factor five. This is possibly a conservative assumption^!, since the 
MECC is expected to have a signal-to-background ratio for this signal slightly better 
than the ECC |117] . We consider, however, this to be the only reasonable choice that 
we can take at this preliminary stage of four-family neutrino detailed study at the 
Neutrino Factory (not to be compared with the order-of-magnitude estimations made 
in previous papers). 

Also in this case, we have grouped events into 10 bins with AEi, = 5 GeV constant 
energy resolution. We have assumed 10% for the bin-to-bin uncorrelated systematic 
errors, 1% and 5% for normalization and energy spectrum distorsion as the correlated 



Notice that, for the silver channel, at least the background induced by right-sign r's with 
wrong charge assignment should be depleted in a four-family neutrino scenario with respect 
to the standard three-family one. This particular background is strongly affected by active- 
sterile mixing angles, since, in the allowed region of the parameter space, — > i^r oscillations 
are significantly depleted with respect to the standard three-neutrino ones. 
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Table 1 

Event rates for the — > Ur and i>e channels for 1 kton MECC detector, exposed to a 

2 X lO^'' (v^, Ve) flux for one year, for different values of On, 624 and 6^4 in the (3+1 ) scheme. 
The other unknown angle, ^13 has been fixed to: ^13 = 5°, 10°. The CP-violating phases are: 
61 = 62 = 0; ^3 = 90°. As a reference, rates at the 1.8 kton OPERA detector (exposed to the 
nominal CNGS beam intensity) and the expected event rates for 1 kton MECC detector in 
the case of the three-family model (i.e., for On = and maximal CP-violating phase 5) are 
also shown. In all cases, perfect efficiency is assumed. 

systematic errors throughout the numerical simulations for both the z/g and the 

^11 signals. Once again, the results here may be slightly conservative because we 

may be overcounting the systematic errors. 

In Tab. [T] we show the expected number of r~ from Vj. and r"*" from z/g 

for a 1 kton MECC detector with perfect efficiency, exposed to a 2 x 10^*^ (z/^, 
flux for one year, for different values of ^13, ^14, ^24 and 6'34. The other parameters 
are: Qx2 = 34°; ^23 = 45°;Am2^i = 7.9 x 10"^ eV^; Aml,^ = 2.4 x lO'^ eV^ and 
AmggL = 1 eV^ (all mass differences are taken to be positive). Eventually, phases have 
been flxed to: 61 = 62 = 0; ^3 = 90°. For comparison, the rates at the CNGS (for the 
nominal CNGS flux, of 4.5 x 10^^ pot/year, an active lead target mass of 1.8 kton 
and 5 years of data taking) and the expected number of events in the three-family 
model for a 1 kton MECC detector with perfect efficiency are also shown. We can 
see that the number of expected r~ events at the 1 Kton MECC is 0(500) at both 
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baselines, with some dependence on the different mixing angles. The fact that at both 
baselines we expect a similar number of events is a consequence of the convolution of 
the — >■ oscillation probability with the z/^iV cross-section and the neutrino 
flux: at the shortest baseline, the probability is maximal below 10 GeV; at the longest 
baseline, the maximum is located in the 30 GeV bin. The higher cross-section for this 
energy bin compensates for the decrease in the i/^ neutrino flux, thus giving a similar 
number of r's in the detector. 

Notice that we will not use the MECC section of the Hybrid-MIND detector to look 
for the golden z/g and the z/^ — >■ z/^ disappearance signals. Both channels can be 

studied, though, at the price of an increase in the scanning load. 



4 Sensitivity to (3 + 1) sterile neutrinos at the Neutrino Factory 

In this section we study the physics reach of the 50 GeV Neutrino Factory setup 
discussed in the previous section to (3+l)-model sterile neutrinos. We will make use 
of all possible oscillation channels available at this setup, namely the golden channel 
T^e ~^ and the disappearance channel z/^ — > using the MIND section of the 
detector; the silver channel Vf. — >■ Uj. and the "discovery channel" z/^ — >■ using the 
MECC section of the detector. 

In Sec. l4.1l and l4.2l we study how a negative result at the Neutrino Factory can be used to 
improve present bounds on the four-family model 9^^'^'^^ , extending the standard three- 
family ^13-sensitivity analysis, and on the three active-sterile mixing angles ^14, ^24 and 
6*34. In particular, the golden and silver channels {ve — ^ t'/x, t-'t) cooperate to put a 
strong bound on ^^3^^™-* (see Sec. 14.11) . On the other hand, the z/^ disappearance channel 
and the "discovery channel" — > v-r improve signiflcantly the exclusion bounds for 
(3-|-l)-sterile neutrinos in the (6'24, 6'34) plane (see Sec. 14.21) . Preliminary results for the 
combination of these two channels at a 50 GeV Neutrino Factory setup have been 
presented in Refs. [TMTTU] . 

In Sec. 14.31 we consider a different situation: in the case of a positive signal at the 
Neutrino Factory, can we discriminate between the (3-|-l)-sterile neutrino model and 
the "standard" three-family one? We will show that, using the combination of the four 
available channels, we can indeed discriminate between the two models in a signiflcant 
region of the active-sterile mixing angles parameter space. In Sees. 14.11 14.21 and 14.31 
we will also compare our results with those that could be obtained using the 20 GeV 
ISS-inspired Neutrino Factory deflned in Sec. El 

In Sec. 14.41 we discuss the dependence of the four channels on the systematic errors. It 
is known that the golden and silver channels are dominated by the statistical errors. 
We will show, on the other hand, that both the disappearance and discovery channels 
are dominated by the systematic errors, emphasizing that improvement on systematic 
errors is extremely important to take full advantage of the discovery channel. 
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In case a positive signal is found, it is interesting to study the various CP-violating 
signals. This is done in Sec. 14.51 where we show the attainable precision in the si- 
multaneous measurement of ^34 and of the CP-violating phase ^3 and the 53-discovery 
potential, when both active-sterile mixing angles 624 and 634 are not exceedingly small. 
We also show how the measurement of the three-family-like CP-violating phase 62 is 
modified by the presence of non-vanishing active-sterile mixing angles. 



4-1 Sensitivity to {9is,9i4) 



In this subsection first we discuss sensitivity to (^13, 614) at the Neutrino Factory. As 
we will see below, sensitivity to 6^14 turns out to be poor. So we will consider sensitivity 
to other combinations of the mixing matrix elements later. 



4.1.1 Sensitivity to {613,614) 

The sensitivity is defined as follows: we first compute the expected number of events 
for z/g —>■ and Ve — >■ Vt oscillations for the input values 6*13^'^™'* = 0, and ^14 = 
^24 = 63A = 0, where 6^^""""^ {{j,k) = (1, 2), (1, 3), (2, 3)) and ^4 = ^-r""^ stand for 
the mixing angles in the four-family scheme denoted by (4fam)lfJ This number, that is 
identical in the three- and four-family models, is labeled as A^". We then compute the 
expected number of events in the (^(3^^™'', 6'i4)-plane for the same oscillation channels 
in the four-family model. The Ax^, computed with respect to the "true" value A^*^, is 
eventually evaluated. The contour for which the 2 d.o.f.'s Ax^ is Ax^ = 4.61 defines, 
then, the region of the parameter space of the (3+l)-sterile neutrino model that is 
non-compatible at 90% CL with the input data corresponding to vanishing {6[^'^^\ 614) 
(to the right of the contour line) and the region that it is still allowed (to the left of 
the line) at this CL. 

The minimum of the is computed as follows: 



Ay = min 

marg par 



— ((1 + a, + a;j/3,)A^j(4fam) 

pol.,(chan.),(L) " '''^ [ j ^3 



+ (1 + ttb + x,/36)S,(4fam) - - Bff + f^- 




+ AXatm+rc(4fam) 



0", 



(9) 



^ Using the four-family expressions for the oscillation probabilities in vacuum shown in 
Ref. [60], it can be seen that, for vanishing 614 and 624, the four- family mixing angle 9^^^^™^ 
maps into the three- family one, 9'f^^^\ On the other hand, for non- vanishing ^14, the four- 
family parameter is expected to be slightly smaller than the three- family one. 
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where we have introduced the prior that comes from the four-family analysis of the 
atmospheric and reactor datalZ 

, , (4^^^-) - 0.50)^ , (|Amy^^-V2.4xlO-3eV^)^ 

AXatn.+.e(4tam; - ^^^^^^^ + ,2(|A^n4fam)|) 



2 (4fam) _ n nn2 (5^4)^ ("^D^ ''"^ ''^ 



'13 - 0.01)" l^i>i4; 1^^,24; , l'^34; 



where | Am32*"^^^™'' | stands for the atmospheric mass squared difference in the four-family 
scheme, and the errors of the oscillation parameters in the four-flavor scheme in eq. 
(fTUj) are deduced from Refs. [21] and [20] as follows: 

a{s\^) = 0.07, ailAml^l) = 0.12 x lO^W, aisj^) = 0.016, 

(j(s?4) = 0.013, a(s^4)=0.02, a(s^4) = 0.12. (11) 

In the minimization procedure in eq. "marg par" stands for the oscillation param- 
eters to be marginalized over (that can be different for different plots), and as, ab, 
(3s and (3b are the variables for the correlated systematic errors, which stand for the 
uncertainties in the overall normalization and in the linear distortion in the spectral 
shape in the magnitude of signal (s) or background (b) [121] . where we have defined 
Xj = Ej/{E^g^^ — -Emin) for neutrino energy Ej for the j-th bin. Following Ref. |121] . 
we assume the correlated systematic errors cTq^ = ctqj^ = 0.01 for the normalization and 
= o"/3b = 0.05 for the spectrum distorsion. In the analysis of the case with single 
baseline length, we have minimized the for each baseline separately, i.e., no sum is 
performed over L, and in the analysis combining the two baselines, we have minimized 
the sum of for each baseline. Similarly, in the analysis of a single channel, no sum 
is performed over the channels ("chan."), while in the analysis combining the different 
channels we have summed over the different channels, i.e., golden and silver in this 
Section. In all cases we sum up for the two possible stored muon polarities ("pol."). 
The index j runs over 10 energy bins. Bj is the background correspondent to the j-th 
bin (5° stands for the expected background in the four-family model for vanishing 
^(4fam) bin). Within this procedure, for any minimization that we 

perform the best-fit for the variables on which we marginalize over can be different. 
However, when we project onto the two-dimensional (^is^'^™^ 6'i4)-plane (or onto the 
(^24; ^34) -plane in the next section) this information is lost. Notice that this approach 
is analogous to that used in the three-flavor analysis in Ref. |120j . A final comment is 
in order: in the case of a real experiment, the minimum of the distribution will in 
general not correspond exactly to the "true" values of the parameters to be fitted. The 



^ Since Ax^ is expected to depend little on the solar neutrino oscillation parameters, we 
will not vary the solar oscillation parameters throughout this paper, so we omit the terms on 
the solar parameters here. 
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Ax^, computed with respect to the best-fit point, will be therefore slightly different 
from our definition. This procedure, however, is commonly adopted to compute the 
expected sensitivity of an experiment that is still under development, so as to compare 
different setups on equal footing. 

The variance is defined as: 



aJ = iV° + 5; + [/,iV0]2 + [/,.5jf , (12) 



where fj is the uncorrelated bin-to-bin systematic error in the j-th bin: 2% for the 
golden channel and 10% for the silver channel, irrespectively of the energy bin, of 
the baseline and of the stored muon polarity. Notice that we can use the GauBian 
expression for the throughout our numerical simulation, with the possible exception 
of the silver channel data at the longest baseline, for which a Poissonian expression 
could be more appropriate due to the extremely low statistics (as it can be seen from 
Tab. [1]). However, as it will be shown in the following discussion on the results, the 
impact of the silver channel at that baseline is negligible. 

In the numerical analysis of this section, the following parameters in common between 



three- and four-family models have been kept fixed to their central values]^ 612 = 34°, 
023 = 45°; Am^^i = 7.9 x 10"^ eV^ Am^^ = 2.4 x lO^^ eV^ The following two 
parameters specific to the four-family model have been also kept fixed: Am|gL = 1 eV^ 
and Si = 0. Eventually, we have marginalized over 6*24 G [0,12°], 9u G [0,35°] and 
^2,^3 G [0,360°]. Matter effects have been included considering a constant matter 
density p = 3.4 g/cm^ for the shortest baseline and p = 4.3 g/cm^ for the longest one, 
computed averaging over the density profile in the PREM [HS] along the neutrino path. 
We have checked that marginalization over a 10% matter density uncertainty does not 
modify our results. 

It is useful to show here the analytic expressions for the golden and silver channels 
oscillation probabilities in vacuum. For both channels it is mandatory to go beyond 
the O(e^) of the expansion in power of small parameters that was shown in eqs. (jUlH]). At 
high orders in powers of e, it is no longer possible to neglect terms proportional to the 
solar mass difference. In particular, for the short baseline L = 3000 km (the only case in 
which vacuum formulae give some insight on the numerical results), the parameter A21L 
ranges from O(e^) to O(e^), depending on the neutrino energy. We will thus show the 
oscillation probabilities Pg/^ and Pgr to order in powers of ^/9l3, v^^, V^^23, 6*34 
and (A2ilv)^/^. For completeness, we have also computed the oscillation probability Pgg 
and Pee to the same order in e {Pes can be found in the Appendix), checking unitarity 
of the mixing matrix. Oscillation probabilities in matter at this order in e are difficult 
to obtain and they will not be presented here. 

The golden channel oscillation probability expanded to order e*^ in vacuum is: 



We will not introduce any label to distinguish these parameters between three- and four- 
family models. 
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+ 2 V2 013^14^24(1 + 5^23) sin (62 -63 + sin 



+ sin 26112 ^13(^21-^) cos [61-52 + 63- ^^^^ sin ^7^^ 



+ ^ sin 2^12 0i4024(A2ii:) sin^i - sl^ el^{A2iL) sin A31L + ^ sin^ 20i2(A2iL)2 , 
V 2 2 

(13) 

whereas the silver channel oscillation probability at the same order is: 



Per — 2 0^,9 2, 

+ 2 {eUl - 26023 - Ol, - el - el + 25^23^34) - ^13^24^34 cos 53} sin^ ^ 

- 2V2 0i30i4024^3'4 Sin (62 -63 + ^) Sin ^ 

A3iL\ . A31L 
sm 



+ 2^2 6^136^14^^34 (1 - 6e23 - ^-f^ sin 1^62 



rN„„„6x X ,X ^31^\ .„ A31L 



- Sin 2012 ^13(1 - di4)i^2iL) COS {6,-62 + 63- sin 
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- sin 2012 0i4034(A2ii:) sin(5i + ^3) - 4 ^^3(^21^) sin A31L + ^ sin' 20i2(A2i 



Notice that the golden and silver channel expressions have leading terms proportional to 
013, i.e. they are 0{e^) (as it was anticipated in Sec. l2.4p . This is the same order at which 
leading terms arise in the three-family model. We expect, thus, a similar sensitivity to 

013 in both three- and four-model analyses at the Neutrino Factory using these channels. 
On the other hand, we expect poor sensitivity to 0i4 using these channels. First of all, 
dependence on 0i4 in the golden and silver channel arise only at higher orders in e (at 
O(e^) for the golden channel and at O(e^) for the silver channel, respectively). Secondly, 
all the terms proportional to 0i4 in the above probabilities appear in combination with 
either 034 or 024, so when we evaluate the value of Ax' by marginalizing over 034 and 
024, both of which have tendency to be small because of the priors, these two probabihty 
can become almost independent of 0i4. For this reason, we expect poor sensitivity to 

014 using these channels. Nevertheless we discuss the sensitivities to both 0i3 and 0i4 
because the sensitivity to 0i3 can have nontrivial dependence on the true value of 0i4 
in the four neutrino case, as we will see below. 

Eventually, the z/g disappearance probability is: 
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Fig. 3. Sensitivity limit at 90% CL in the (sin^ 2^13, 614) plane for 63 = and different values 
of 62- The solid lines refer to the golden channel results, only. Dashed lines stand for the sum 
of golden and silver channel results. The colors are: blue for L = 3000 km; red for L = 7500 
km; green for the combination of the two baselines; the horizontal dashed grey line represents 
the present bound on 614. The four panels represent our results for: 62 = (top left); 82 = 90° 
(top right); 82 = 180°; (bottom left); 82 = 270° (bottom right). 

Pee = 1 - 2 eUl - el) - AoUi - el, - 2 el) sin^ ^ 

+ 2 sleK/\2iL) sin A31L - sin^ 2^i2(A2iL)2 (15) 

Notice that this channel has leading O(e^) dependence on both ^13 and e^. For this rea- 
son the most stringent bound on this angle has been obtained by reactor experiments, 
as we have shown in Sec. 12. 3[ A detector with an extremely good electron identification 
efficiency is needed to perform this task, however, something beyond the reach of the 
MIND or MECC detectors. We are, thus, not considering this channel in the present 
paper. 

Figs. [3] and m present the sensitivity to (3+l)-sterile neutrinos in the {e^ii^^\ 6'i4)-plane 
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Fig. 4. Sensitivity limit at 90% CL in the (8111^2013, 6*14) plane for 83 = 90° and different 
values of 62- The solid lines refer to the golden channel results, only. Dashed lines stand for 
the sum of golden and silver channel results. The colors are: blue for L = 3000 km; red for 
L = 7500 km; green for the combination of the two baselines; the horizontal dashed grey line 
represents the present bound on ^14. The four panels represent our results for: 82 = (top 
left); 62 = 90° (top right); 62 = 180°; (bottom left); 62 = 270° (bottom right). 

at 90% CL at the 50 GeV Neutrino Factory setup for several representative choices of 
62 and 63. In particular, ^3 has been fixed to ^3 = and 90° in Fig. [H] and Fig. HI 
respectively. In both cases, the phase 62 (that reduces to the three-family CP-violating 
phase 6 in the limit da — >■ 0) has been fixed to 62 = 0, 90°, 180° and 270°. Results 
for the two baselines are shown both separately and summed: blue lines stand for the 
shortest baseline data; red lines stand for the longest baseline data; green lines stand for 
the sum of the two baselines. Solid lines stand for golden channel data, only, whereas 
dashed lines stand for the sum of golden and silver channels data. Eventually, the 
horizontal dashed grey line represents the present bound on 9u that takes into account 
the four-family analysis of the atmospheric and reactor data. 

First of all, notice that the golden and silver channel data show a very limited sensitivity 
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to 9i4 when marginalizing over and ^34, as it was expected by inspection of the 
oscillation probabilities Pe^ and Per, eqs. fll3fl4p . 

Secondly, when we focus on the golden channel results (solid lines), we notice a rather 
different behavior at the short and long baseline. At the short baseline, the sensitivity 
to 61^^^^ is significantly phase-dependent: the maximal sensitivity to ^^3^^™'' ranges 
from sin^ 29^^^'"^^ = 1.5 x IQ-^ for ^2 = ^3 = to 6 x IQ-^ for 52 = ^3 = 90°, with 
a strong dependence on the 6*14 value. This is a consequence of cancellations between 
the term proportional to 6'i3(A2i-L) in the third line of eq. ffT^ . that in three-family 
represents the interference between solar and atmospheric oscillations, and the term 
proportional to 6'i36'i4^24 in the second line of the same equation. 

The behavior is extremely different when we consider the golden channel at the long 
baseline. The location of the far detector corresponds to the three-family "magic base- 
line" |105j . where the three-family dependence on the CP- violating phase 6 vanishes. 
Notice that this happens as a consequence of the vanishing of the term in the third 
line of eq. ffT^ when computed in the Earth matter. This is, indeed, the only term 
through which a ^-dependence enters into the three-family golden channel oscillation 
probability. The four-family term in the third line of eq. ( fT3|) . that reduces to the 
three-family one for vanishing ^,4 and 61,63, vanishes also. Therefore, no cancellations 
between different terms occur at this baseline and sensitivity to 9{^^^^ depends much 
less than for the short baseline on the values of 62 , 63 and ^14 . This is precisely what 
we can see in all panels of Fig. [3] and HI 

A similar effect can be observed when we add the silver channel data to the golden 
channel ones at the short baseline. In this case, the term in the third line of eq. ( fT3l) 
cancels with the term in the fifth line of eq. ( IT^ at the leading order, leaving an 
O(e^) term suppressed by 9r39l^{A2iL). Also in this case, we see indeed in all panels 
of both figures that a reduced dependence of the sensitivity to 6^(3^'^™^ from 62, 63 and 
^14 is achieved. Particularly striking is the case of 62 = 63 = 90°, Fig. H] (top right), 
where we can see that the sensitivity to ^^3^'^™'* goes from sin^ 29[^^^^ = 6 x 10~^ for the 
golden channel alone to sin^ 29^^^^^ = 3 x 10~^ for the combination of golden and silver 
channels. On the other hand, the silver channel statistics at L = 7500 km from the 
source is extremely poor (see Tab. [1]). For this reason we have no impact whatsoever 
of this channel on the 9[^^^^ sensitivity at the long baseline. 

The performance of the detectors located at the two baselines differ significantly de- 
pending on the specific choices of 62 and ^3. For most of the choices of (^2, ^3), the longest 
baseline outperforms the shortest one, with the only exceptions of 62 = 0, 180°; ^3 = 
(Fig. El left panels) and 62 = 180°; ^3 = 90° (Fig. H bottom left panel). 

The combination of the two channels and the two baselines at the 50 GeV Neutrino 
Factory has a rather good sensitivity to 9[^3^^^: in the absence of a signal, we can 
exclude values of ^^3^*^™^ larger than sin^ 29[^^^^ < 4 x 10~^ for any of the choices of 
{62, 63) shown. The sensitivity to the active-sterile mixing angle ^14, on the other hand, 
is 9i4 < 9° when marginalizing over 6'24 and ^34. 
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Fig. 5. Sensitivity limit at 90% CL in the (sin 2^13, 614^) plane, marginalizing over 624,, ^34, ^2 
and 83. The solid lines refer to the golden channel results, only. Dashed lines stand for the 
sum of golden and silver channel results. The colors are: blue for the shortest baseline; red for 
longest baseline; green for the combination of the two baselines; the horizontal dashed grey 
line represents the present bound on 614. Left panel: 50 GeV Neutrino Factory; Right panel: 
20 GeV Neutrino Factory. 

We eventually compare the results obtained using the 50 GeV Neutrino Factory setup, 
Fig. Et^left), with those that can be achieved using the 20 GeV Neutrino Factory ISS- 
inspired setup, Fig. [5](right). As before, solid lines stand for golden channel data, only, 
whereas dashed lines stand for the combination of golden and silver channel data. Blue 
lines stand for the shortest baseline (3000 km or 4000 km); red lines for the longest 
baseline (7500 km in both setups); green lines for the combination of the two; eventually, 
the grey dashed line represents the present bound on Ou. In addition to 6*24 and 6*34, we 
have marginalized over the CP- violating phases 62 and ^3, with 62,63 G [0,360°]. The 
rest of the parameters have been kept fixed to the values given above. 

First of all notice that for both setups the ^j3^^™''-sensitivity that can be reached using 
golden channel data, only, is much poorer at the short baseline than at the long baseline. 
We have found O^^^^^^ < 6 x 10~'*(9 x 10~^) for any value of 614 at the short baseline 
for the 50 GeV (20 GeV) setup, to be compared with ^Jg^^"^^ < 2 x 10-^(3 x 10"^) 
at the long baseline for the 50 GeV (20 GeV) setup. The second relevant point is 
that the silver channel significantly improves the ^[g^'^^'^-sensitivity for the 50 GeV 
setup, but only marginally for the 20 GeV one. We can see that the ^|3^'*™''-sensitivity 
at the short baseline using the combination of golden and silver channels becomes 
^(4fam) ^ 2 X 10-^(6 X 10"^) for the 50 GeV (20 GeV) setup. When considering the long 
baseline we can see, however, that the impact of the silver channel becomes marginal 
(though it can still be observed for the 50 GeV setup). 

The ^[g^'^^'^-sensitivity that can be reached using the combination of the two baselines 
is: ^S^"^^ < 6 X 10-5(1.5 x 10"^) for any value of 9u at the 50 GeV (20 GeV) setup. 
The inclusion of the silver channel data does not modify significantly these results. 
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Our conclusion is that the higher energy setup has a much greater ultimate sensitivity 
to eif''"'^ than the ISS-inspired 20 GeV one. On the other hand, both setups are not able 
to improve the present bounds on 6*14 when marginalizing over ^24,^34 (see, however. 
Sec. 14. 3p . The silver channel can significantly improve the 6'i3^^™''-sensitivity when only 
one baseline is considered at the 50 GeV setup. However, it has a negligible impact 
when the combination of the golden channel data at the two baselines is considered. 

Eventually, we have studied the impact of the correlated systematic errors as, /3s in 
eq. ([9]) on our results. For both the 50 GeV and 20 GeV setups, we have found that 
the results shown in Fig. [3] do not change for = /Sg = 0. Golden and silver channels 
are indeed dominated by statistical errors. 



4- 1-2 Sensitivity to Ue^U^^ and UeiUrA 

Since sensitivity to 614 at the Neutrino Factory is poor, it is worth investigating whether 
the Neutrino Factory has sensitivity to other combinations of the mixing matrix ele- 
ments. From the form of the appearance oscillation probability 



P{i^a ^ 1^(3) = 4Re U^3U;,{U:,Ups + K.U^,] 



sm^ ' 



AE 



we can expect that the golden and silver channels have some sensitivity to UeAUn^ and 
UeAUrA- In the prcscut parametrization ([2]) of the mixing matrix, we have Ue^Un^ = 
S14C14S24 = S14S24 + 0{e^) and Ue^Uri = S14C14C24S34 = S14S34 + O(e^), where we have 
used the fact that \9i4\ and 16*241 are small. In the analysis of sensitivity to f/e4f/^4 
{Ue4Uri), instead of using \9i4\ and |^24| (I ^14 1 and 16*34!) as the independent variables, 
it is convenient to take S14S24 and S14/S24 (S14S34 and S14/S34) as the independent ones, 
respectively. Taking this basis we have performed analysis on sensitivity to [/e4f/^4 and 
UeiUri- The sensitivity is defined as in Sect. 14.1.11 with straightforward replacements. 
The results are shown in Fig IH] and they indicate that the combination of the golden and 
silver channels has good sensitivity to these variables. To see how much improvement 
we have, we take the square root of t/e4?7^4 or UeiUri, so that these factors correspond 
roughly to sine of some angle. Then the upper bound for ^JUe/JJ^^ by the 50 GeV (20 
GeV) Neutrino Factory is VSxTcR (Vl x 10-3), which is about 15% (20%) of the 
current bound, ^JUeiU^ < VO.02. Similarly, the upper bound for y/Ue4Ur4 by both 50 

GeV and 20 GeV Neutrino Factory is v2xl0^, which is about 20% of the current 
bound, VUIJJ^ < 0.06. 

For both UedJ^^ and Uedlr^ plots, we see that both the golden and silver channels play 
a role in giving the constraints. As it is expected from statistics, the result by the 50 
GeV Neutrino Factory is better for UedJ^^, but the sensitivity to UeiUr^ is almost the 
same for the two setups (notice that the scale on the vertical axis for the left and right 
panels are different). For the 20 GeV case, the data at 7500 km perform very well and 
the combined data of 4000 km and 7500 km give a result almost comparable to that of 
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Fig. 6. Sensitivity limit at 90% CL to UeiU^i and Ue^UrA- Ax^ is evaluated for a fixed set 
of values of (sin^ 2^13, UeiU^i) (upper panels) or a fixed set of values of (sin^ 2^13, UeiUrA) 
(lower panels), marginalizing over S14/S24, 6*34, 82 and ^3 (upper panels), or over S14/S34, 
^24; ^2 and (53 (lower panels). Left panels: 50 GeV Neutrino Factory; Right panels: 20 GeV 
Neutrino Factory. The current bound on UeiU^i (UeiUri) is 0.02 (0.06). 



50 GeV. 

It is interesting to note that the golden channel also plays a role in improving sensitivity 
to UeiUri oc S34. We have obtained a lengthy analytic formula for the golden channel 
in matter with some approximation, i.e., to quartic order in e while keeping all orders 
in ^34 and we have found that dependence on 6*34 appears through the form of 0(e^) x 
s\^Ae/ ^E'ii. This explains why the golden channel has some sensitivity to UeJJri 
through the matter effect at the longer baseline, once we choose a suitable set of the 
independent parameters to vary. 
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4-2 Sensitivity to (^24, ^34) 



The sensitivity is defined as in the previous section: we first compute the expected 
number of events for and oscillations for the input values Ofl'^^^ = 

and 9u = 6*24 = 6*34 = 0, N^; we then compute the expected number of events in the 
(^24, ^34) -plane for the same oscillation channels in the four-family model. The Ax^, 
computed as in eq. with respect to the "true" value N^, is eventually evaluated. 

Contrary to the case of the golden and silver channels, however, the number of expected 
background events is much smaller than the signal for both the disappearance and 
discovery channels. The effect of and jSb is, thus, negligible. Hence, we will not 
perform minimization with respect to and /3f, for these channels and put these 
parameters to zero hereafter. As in the case of the golden and silver channels, the 
variance is defined by eq. ( fT2l) . where Nj in this case are the number of events of 
the disappearance or discovery channels, and the uncorrected bin-to-bin systematic 
error fj is 5% for the disappearance channel and 10% for the discovery channel, 
irrespectively of the energy bin, of the baseline and of the stored muon polarity. 

In the numerical analysis, the central values of the following parameters in common 
between three- and four-family models are: 6*12 = 34°, 6*13 = 5.7°, ^23 = 45°; Am^^i = 
7.9 X 10-5 eV^ Aml^ = 2.4 x 10"^ eV^; ^2 = ^ = (where S is the three-family CP- 
violating phase). The central values for the following three parameters specific to the 
four-family model are: Arngg^ = 1 eV^, 6*14 = and Si = 0. Matter effects have been 
included considering a constant matter density p = 3.4 g/cm^ for the shortest baseline 
and p = 4.3 g/cm^ for the longest one, computed averaging over the density profile in 
the PREM [86j along the neutrino path. We have checked that marginalization over a 
10% matter density uncertainty does not modify our results. 

The sensitivity to (3+l)-sterile neutrinos at the 90% CL in the (6*24, 6^34) -plane for the 
50 GeV Neutrino Factory setup is shown in Figs. [7] and [H In these figures, we have 
studied how the marginalization over ^23,^^31 and ^3 modify our results by varying 
them in the ranges ^23 e [40°, 50°], Amj^ e [2.0,2.8] x lO'^eV^ and ^3 G [0,360°]. 
We have also checked the impact of the marginalization over all the other parameters 
(otherwise fixed to their central values, given above), by studying them one by one in 
combination with 6'23, Am|j^ and ^3. The considered marginalization ranges are: 612 G 
[30°, 36°]; 013 e [0,10°]; 014 G [0,10°];Amii G [7.0,8.3] x 10" W; 5i G [0,360°] and 
S2 G [0,360°]. We have found that none of these parameters modify significantly our 
results, contrary to the case of 6'23,Am3i and ^3. Eventually, we have checked that 
changing the sign of Am|i does not modify our results, either] ^^1 

First of all, in Fig. [7] we show the sensitivity limit at 90% CL to 6*24 and 6*34, for fixed 
63 = 0, whilst marginalizing over 6*23 in the range 023 G [40°, 50°] and Am3i in the 
range Am|^ G [2.0, 2.8] x 10^^ eV^. Notice that the considered allowed range for 023 is 
a bit smaller than its present allowed range from the three-family global analysis. The 

Notice that we are not sensitive to the sign of the largest mass difference, Am|i. 
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Fig. 7. Sensitivity limit at 90% CL to 624 and 634, for fixed 63 = 0, whilst marginalizing over 
023 in the range 623 G [40°, 50°] and Am^^ in the range Am|^ E [2.0,2.8] x lO'^ eV^. Red 
lines stand for the disappearance channel — > i/^; hlue lines stand for the discovery channel 
i/^ — > i/^; green lines stand for the combination of both; the dashed grey line represents the 
present bound on 624 and 634. Left: L = 3000 km baseline, Right: L = 7500 km baseline. 

Neutrino Factory, however, has an enormous potential for improving the measurement 
of the three-family atmospheric parameters through the z/^ disappearance channel as it 
was shown, for example, in Ref. P-09]. It is, therefore, absolutely reasonable to vary 6*23 
over a reduced range. In the two plots, red lines stand for the — >• z/^ disappearance 
channel data; blue lines stand for the Vr discovery channel data; green lines stand 
for the combination of both; the dashed grey hne represents the present bound on 6*24 
and 6*34. 

We discuss first the left panel of Fig.[71 that refers to the results at the shortest baseline, 
L = 3000 km. Notice that the disappearance channel has the strongest sensitivity to 
^24, that can be excluded above ^24 ~ 8° for any value of ^34. On the other hand, 
the Ufj^ Vr channel gives the ultimate sensitivity to ^34 for vanishing ^24, ^34 < 
15°. In both channels, we can see a strong correlation between ^24 and ^34, with a 
bound on 6*34 that is strongly dependent on the specific value of 624 considered. This 
behavior can be understood by looking at eqs. (p[7|) in Sec. [21 The strong correlation 
between the two mixing angles is induced by the subleading O(e^) terms, proportional 
to (A„L)s24S34 cos53 (siucc we are considering ^3 = the term proportional to sm.53 
in P^r vanishes). The (6'24, 6'34)-correlation in the disappearance channel is softer 
than in the z/^ one: a consequence of the different statistical significance of this 

term in appearance and disappearance channels. Eventually, the combination of the 
two channels gives a very good sensitivity to both mixing angles. 

A similar combined sensitivity is achieved at the longest baseline, L = 7500 km, whose 
results are shown in the right panel of Fig. [71 We can see comparing the blue lines 
on both panels that the z/^ —>■ v-r channel sensitivity to (^24, ^34) is substantially the 
same when changing baseline (the expected number of events for this channel at the 
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two baselines is very similar, see Tab. [T]). On the other hand, the z/^ disappearance 
channel data (the red line) has a much stronger sensitivity than at the short baseline: 
a consequence of the increased significance of the subleading term with respect to the 
dominant one in the disappearance channel due to larger matter effects at 7500 km. 

The sensitivity of the disappearance channel to du arises at a higher order in the 
expansion that we have presented in eqs. (pfTl) . If we introduce terms to fourth order 
in e (under the assumption 613 = 6*14 = 0, and taking into account the deviations from 
maximality of 623), we get: 
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Most of the 6'34- dependent terms in eq. (1161) are proportional to the matter parameter 
(AnL). This means that the impact of these terms will be more important at the longest 
baseline, than at the shortest one (as we have seen in Fig. [7|). On the other hand, the 
^24-sensitivity arises from the 6*34 term at 0{e'^) in eq. ffTUj) . 

These behaviors are strongly modified if we marginalize over the CP-violating phase 
^3, as it can be seen in Fig. [HI where we present the sensitivity limit at 90% CL to 
^24 and ^34, whilst marginalizing over ^23 in the range 6*23 G [40°, 50°], over ^3 in the 
range ^3 G [0, 360°] and over Aml^ in the range Am|i G [2.0, 2.8] x lO^^eV^. As before, 
red lines stand for the disappearance channel data; blue lines stand for the 

discovery channel data; green lines stand for the combination of both; the 
grey dashed line represents the present bound on 6^24 and 6*34. 

Notice that the correlation between ^24 and ^34 in the disappearance data (red 
lines) has vanished. This is a straightforward consequence of the marginalization over 
^3, that removes the term proportional to (A„L)s24S34 cos ^3 in eq. (Ej), responsible for 
that correlation. We have found that the sensitivity on the two active-sterile mixing 
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Fig. 8. Sensitivity limit at 90% CL to 824 and 634, marginalizing over 623 G [40°, 50°], the 
CP-violating phase S3 G [0,360°] and Am|i E [2.0,2.8] x lO'^eV^. Red lines stand for the 



disappearance channel Uu 



blue lines stand for the discovery channel — > Vt-; green 



lines stand for the combination of both; the grey dashed line represents the present bound 
on 824 and 834. Left: L = 3000 km baseline, Right: L = 7500 km baseline. 

angles 624 and 6*34 is now represented by vertical and horizontal lines (similar to the 



results in Fig. A similar effect is found for the 



V-r 



appearance channel data 



(blue lines), though softer. We can see that, at both baselines (but more significantly 
at the longest one) some remnants of the (6*24, 6*34 )-correlation can still be observed 
for this channel. Comparing the results of the two channels, we have found again that 



at the short baseline the v.. 



data give the best sensitivity to 6*24, whereas the 



^11 i^T data give the best sensitivity to ^34. At the long baseline. 



z/^ is as good 



as Uj- in the 6'34-sensitivity. 



We, eventually, compare the 90% CL (6*24, 6'34)-sensitivity that can be obtained using 
the combination of the two channels and of the two baselines at the 50 GeV Neutrino 
Factory setup, Fig.[9]^left), with the one that can be achieved using the 20 GeV Neutrino 
Factory ISS-inspired setup, Fig.[9](right). As before, red lines stand for z/^ disappearance 
channel data; blue lines for the z/^ Vj. discovery channel data; green lines for the 
combination of both channels; the grey dashed line represents the present bound on 
624 and ^34. Data for the two baselines are always summed. In these plots, we have 
marginalized over 623 e [40°, 50°], ^3 e [0,360°] and Amii G [2.0,2.8] x IQ-^eV^ The 
rest of the parameters have been kept fixed to the values given previously!^ 

First of all, we can see by comparing red lines (the disappearance channel data) be- 
tween Fig. [9]^left) and Fig. [fright) that the ultimate sensitivities to 624 and 6'34 at the 
two setups are very similar: the upper bounds 624 < 7.5°(8°) and ^34 < 15° can be 
inferred from the data for the 50 GeV (20 GeV) setup. When we eventually combine 
the results for the disappearance and the discovery channels, however, we find that the 



We have checked that the effect of the marginalization on the rest of the parameters do 
not affect the results. 
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Fig. 9. Sensitivity limit at 90% CL to 624 and 634, marginalizing over 623 G [40°, 50°], the 
CP-violating phase 63 £ [0,360°] and Am|i e [2.0,2.8] x 10"^eV^, for the combination of 
the two baselines. Dashed lines stand for the results without the correlated systematic errors. 
Red lines stand for the disappearance channel v^; blue lines stand for the discovery 



channel — > Vr', green lines stand for the combination of both. Left panel: 50 GeV Neutrino 
Factory; Right panel: 20 GeV Neutrino Factory. 



50 GeV Neutrino Factory outperforms the 20 GeV ISS-inspired one, as it can be seen 
comparing the green hues in Fig. [9l This can be easily explained pointing out that the 
discovery channel data (blue lines) are able to exclude a significantly larger region of 
the parameter space when going to higher energy, a straightforward consequence of the 
higher statistics due to the higher Ut-N cross-section. 

In Fig. [HI we have also studied the impact of the correlated systematic errors as, Ps 
in eq. © on our results. Dashed blue, red and green lines represent the 90% CL 
sensitivities to (6^24, 6^34) when the correlated systematic errors as,Ps are not included 
and only the uncorrelated bin-to-bin systematic errors {fj = 5% for the z/^ disappearance 
channel and fj=10% for the discovery channel, irrespectively of the energy bin, of the 
baseline and of the stored muon polarity) are taken into account. 

Comparing dashed and solid lines at the 20 GeV Neutrino Factory (right panel), we 
can see that the inclusion of the correlated systematic errors has a negligible impact 
when analyzing the data at this setup. On the other hand, when studying the data for 
both channels at the 50 GeV setup, we can see that the inclusion of these errors modify 
our results. In particular, the sensitivity to 6*24 through the z/^ disappearance channel 
goes from 6'24 < 6° when only fj is considered to 6^24 < 7.5° when a^, Ps are also taken 
into account. 

Our final conclusion is the following: using the 20 GeV ISS-inspired Neutrino Factory, 
we get two roughly independent limits on the two angles 6*24 and 6*34 : 6*24 < 8° for any 
value of 6'34 and 6'34 < 14° for any value of 6^24 at the 90 % CL. Slightly stronger ultimate 
limits are obtained at the 50 GeV Neutrino Factory: 6'24 < 7.5° for vanishing ^34 and 
^34 < 12° for vanishing ^24- The significantly large discovery channel statistics at the 50 
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GeV Neutrino Factory with respect to the 20 GeV one, however, strongly increase the 
sensitivity of this setup to the combination of 6*24 and 6*34, such that a roughly diagonal 
line in the (6*24, 6*34) plane connecting (^24,0) and (0,6*34) can be drawn. 



4-3 Discrimination of the four neutrino schemes 



In the subsections 14.11 and 14.21 we have discussed the sensitivity to 6*13, 6*14, 6*24, 634 by 
looking at statistical significance of deviation of a four-flavor scheme from that with 
a certain set of reference values of the oscillation parameters. Here we will discuss 
whether the Neutrino Factory setup can distinguish our four-neutrino scheme from the 
three-flavor scenario. 

We introduce the " sterile neutrino discovery potential" , defined as follows: 



AY^(4fam) = min 

marg par 



,pol.,(chan.),(L)"'''^' [ j ^3 

+ (1 + ab + Xjpb)Bj{3iam) - Ar,(4fam) - 5j(4fam))^ 

V I , A 2 




+ ^Xatm+re(4fam) ] , 
(19) 



where Ax^^jj^^j.g(4fam), defined in eq. ( fTOl) . is the prior from the four- flavor oscillation 
analysis of the atmospheric and reactor data, and the errors of the oscillation parame- 
ters in the prior Ax^tjjj_,_j,g(4fam) are given by eq. (fTTi) . 

A remark is in order. The definition of the Ax^ in the present case, although looking 
similar, is slightly different from that used in the previous sections. In Sees. 14. II and [4.21 
we assumed that the minimum of the corresponds to the "true" values of the four- 
family model, and therefore Xmm 4fam — 0- The Ax^ is then computed in the same model 
in which data are generated, and CL contours define the region of parameter space 
compatible at a given CL with the "true" values. In this Section we also assume that 
data are generated in the four-family model, but we try to fit them in the three-family 
model. The minimum of the in the four- family model is located at the "true" values 
of the parameters and Xmin 4fam — 0- On the other hand, when we try to fit the four- 
family-generated data in the three-family model we will in general find Xminsfam 0; 
since a "wrong" model is used to fit the data, except for the special case defined by 
9,, = (j = 1,2,3), C^"^) = %|bestfit ((^,j) = (1,2), (1,3), (2, 3)), Am'^,^''^^^ = 
^^jilbestfit (j = 2,3), where the two models coincide and A^^ = 0. In the rest of 
the four-flavor parameter space, the Ax^ defined in eq. (|T^ corresponds to Xmin sfam ~ 
Xmin4fam- C!L coutours define, then, regions in the four-family parameter space for 
which a three-family fit to the data is worse than a four-family fit to the data of a 
quantity Ax^. For example, a point with Ax^ = 4.61 is a point that is fitted by the 
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four-family model much better than by the three-family model. We will define points 
outside this contour as points for which the hypothesis that data can be fitted in the 
three-family model is "excluded at 90% CL". Under these premises, we can use eq. ( IT^ 
to determine regions in which we are able to distinguish four- from three-family models 
in the four-flavor parameter space in the same manner as in the previous subsections. 

Since the excluded region is expected to depend little on the solar neutrino oscillation 
parameters in the three-flavor scheme, we will not marginalize Ax^ with respect to 
^(3fam) Am2i*'^^'^™'* in eq. fll9p . Moreover, while Ax^ is naively expected to depend 
on all of the parameters 6^i^'^^\ 9^^^^\ \Am1i^^^^™^\, we have found numerically 

that it suffices to vary some of the parameters and put other parameters to the best-fit 
values in most analyses. Namely, in the case of the golden and silver (disappearance 
and discovery) channels, a dominant role is played by 6^^^^^^ and ^(^f'^™) (^'23^'^™^ and 
I Am3i*'^^^'^''|) which are the only three-family parameters that we vary, and the other 
three-family ones are fixed in the analysis. 

To compare the results with those in the subsections 14.11 and 14.21 we project the ex- 
cluded region either in the {e[f'^\eu) plane, or in the {624, Gsa) plane. In these projec- 
tions, we would like to obtain the most conservative excluded region, i.e., the common 
excluded region in the (6^24, 6^34) plane irrespective of the values of 6^(3^'^™'' and 6^14, or 
the common excluded region in the {0i^^^\6i4) plane irrespective of the values of 624 
and 6*34. Notice that the four angles can in principle be measured simultaneously if we 
use informations from the four channels at the same time. To obtain them, in princi- 
ple we have to marginalize Ax^ not only with respect to the three-family parameters 
described above but also with respect to the four-family ones, such as 0[^'^^\6^^^\ 
Amli ^^^"'\ |Am32^^^'^™''|, 61, 62, S3, as well as (6*24, ^34) in the former case and {Oi^^'^\ ^14) 
in the latter. In marginalizing over the four-family parameters, however, we do not have 
to vary all the parameters for a couple of reasons. First of all, since the excluded re- 
gion is expected to depend little on the solar neutrino oscillation parameters in the 
four-fiavor scheme, we can fix the solar parameters 6[^'^^\ Am^i ^^ , Si. Secondly, be- 
cause of the prior Axit^+j.ei^fa.m) , in practice we can fix the following parameters 

to the best fit values: s^^''""^ ^ 0.01, 5^^"'""^ ^ 0.5, |Am^/''"'")| ~ 2.4 x lO'^eV^, 
(^24, ^34) ^ (0,0) in the case of the (^if """"^ ^14) plane, and (^ir°'\^i4) ^ (5.7°, 0) in 
the case of the {624, O34) plane. Thus, the only non-trivial four-family parameters to be 
marginalized over are S2 and ^3 in the case of the {0^^^^\ 614) plane, and ^3 in the case 
of the (6*24, 6*34) plane. 

The results obtained are presented in Fig. [101 where the dashed black line stands for the 
region which is excluded by the prior, i.e., by the present data of the atmospheric and 
reactor experiments. Upper panels show the "sterile neutrinos discovery potential" of 
golden and silver channels; lower panels the "discovery potential" of disappearance 
and z/^ Ut appearance channel. On the left, we show results obtained for the 50 GeV 
Neutrino Factory; on the right, using the 20 GeV ISS-inspired setup. We can see in 
Fig. fTOlf upper panels) that at both setups, if we have no information on 6'24 and ^34, 
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Fig. 10. Left-upper panel: The right upper part of each Une is the region projected onto the 
(Gi3, 0u) plane, in which the hypothesis of the three flavor scheme is excluded at 90% CL 
at the 50 GeV Neutrino Factory. It is obtained by marginalizing over 6^r^^^^ and ij^^^'^™) as 
well as 62 and 63. The grey dashed lines stand for the excluded region obtained only from 
the prior AXatm+re (^fam) where terms other than (s^g*-'*^'*™-' - 0.01)^/cr^(sf3) + (sl^)"^ / a"^ (sl^) 
are assumed to be zero in eq. ilO\) . The excluded regions for nonvanishing 62A or 634, which 
are always larger than the case for 62A = ^34 = 0, are also depicted for information. Left- 
-lower panel: Excluded region at 90% CL projected onto the (624,, Osa) plane. It is obtained 
by marginalizing over 9^^'^^\ | Amg;^^^^^™^ | as well as 63- The grey dashed lines stand for 
the excluded region obtained only from the prior Axatm+rel'^fam) where terms other than 

S34) are assumed to be zero in eq. ^U^. Right-upper (lower) panel: 
The same hgure as the left-upper (lower) panel for the 20 GeV ISS-inspired setup. 



then the golden and silver channels cannot discriminate between the three- and four- 
family models for both setups within the presently allowed region. If we can previously 
measure 6*24 and 6*34 using the z/^ and channels, finding that one or both 

angles are non-vanishing, then we can see in both upper panels of Fig. [10] that the 
" sterile neutrino discovery potential" of the combination of golden and silver channels 
is significantly improved. For example, for 6'34 ~ 10°, the 50 GeV (20 GeV) setup can 
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distinguish the (3+l)-neutrino model from standard three-family oscillations for On as 
small as ^14 ~ 2° (4°). This is because of the same reason that the Neutrino Factory 
has good sensitivity to UeiUfj^^ and UeiUrA, as was discussed in sect. I4.1.2[ We point 
out that the silver channel gives a non-trivial contribution to the " discovery potential" 
Ax^ in some regions of the parameter space, although the contributions of the golden 
and silver channels are not shown separately. 

In Fig. [TOTlower panels) we can see that, irrespectively of prior knowledge of 6[f^™^ 
and 614, the combination of the z/^ disappearance and the z/^ — > Uj. discovery channels 
permits discrimination between the four- and three-neutrino oscillation models in a 
significant region of the presently allowed parameter space. The additional information 
from golden and silver channels increases the region in which discrimination is possible. 

We conclude that combination of the four channels is extremely effective to tell the 
difference between the four- and three-flavor schemes in a significant region of the 
presently allowed parameter space. Notice that the synergy between the four chan- 
nels is not symmetric: whereas a previous knowledge of (6*24, 6*34) strongly increase the 
"discovery potential" in the {9i^^™\6i4)-p\a.ne, the measurement of (^'13^'^™'*, ^14) has a 
much smaller impact in the "discovery potential" in the (6'24, 6'34)-plane. Eventually, 
the 50 GeV setup has a greater "sterile neutrino discovery potential" than the 20 GeV 
one. 



4-4 Dependence of sensitivity on the systematic errors 

We have also investigated the dependence of the performance of the four channels on 
the systematic errors. 

As for the golden and silver channels, for which statistical errors are dominant, we 
have found from numerical calculations that the sensitivities of these two channels to 
61^^^^ and 614 depend only to some extent on cJo (the correlated systematic error on the 
overall normalization), and they depend very little on fj (the bin-to-bin uncorrelated 
error) and ap (the correlated systematic error in the linear distortion of the spectral 
shape) . 

The dependence of the discovery and disappearance channels on fj, and ct/j (and 
on the MECC volume, in the case of the former) is shown in Fig. [TT] The upper 
panels present the discovery channel (6'24, 6'34)-sensitivity, where we consider fj = /^^ = 
0.10, 0.03, aa = 0.100, 0.050, 0.025, 0.010, and the MECC mass either 4 or 8 kton. The 
lower panels present the disappearance channel (^24, ^34) -sensitivity, assuming fj = 
f^^ = 0.05,0.02 and = 0.100,0.050,0.025,0.010. We have checked numerically that 
the impact of the systematic error in the linear distortion of the spectral shape 
is small and it will not be discussed here. For both channels we have considered the 
50 GeV setup performance (left panels) and the 20 GeV setup one (right panels). 
From Fig. [TT]( upper panels), we see that /^^ is the most important factor to improve 
the performance of the discovery channel for both the 50 GeV and 20 GeV neutrino 
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Fig. 11. The 90%CL sensitivity to (624, ^34) using the i/^ v^- (the upper panels) and — > 
(the lower panels) channels for different values of the uncorrelated hin-to-bin systematic error 
fj = /^T- = 0.1, 0.03 and fj = ff^^ = 0.05, 0.02, of the correlated systematic error on the overall 
normalization = 0.1,0.05,0.025,0.01 and of the MECC mass (= 4, 8 kton in the case 
of Vr)- Left panels: at the 50 GeV setup; Right panels: at the 20 GeV ISS-inspired 

setup. In all the hgures black lines stand for the excluded region for the reference values 
used in the calculations in other sections. The grey lines stand for the excluded region using 
the disappearance (discovery) channel with f^^ = 0.05, = 0.01 (with 4 kton, Z^,- = 0.1, 

CTa, = O.Olj. 

factories. On the other hand, an increase of the MECC mass from 4 to 8 kton improves 
only marginally the discovery channel sensitivity. In Fig. ITlT lower panels), we see that a 
reduction of both /^^ and aa should be pursued to increase the disappearance channel 
sensitivity. One important conclusion from Fig. [H] is that an improvement of f^r 
below 10% is mandatory in order to take full advantage of the discovery channel at 
the Neutrino Factory, particularly for the 20 GeV setup. This error represents indeed 
our incomplete knowledge of the MECC detector. To improve our analysis on the 
discovery channel in the future, we need detailed information on the correlated and 
uncorrelated systematic errors, such as the uncertainties of the detection efficiency 
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which depend only on the nature of MECC or depend on the characteristics of the 
individual detectors located at each of the two baselines. Although there has been no 
study on these systematic errors so far, they are expected to be better understood 
after the first years of data taking of the OPERA experiment (that started operation 
in 2008). 



4-5 A CP-violating sterile neutrino signal 



In Sees. 14.11 and 14.21 we have considered the case of a null result for sterile neutrino 
searches at the Neutrino Factory after 4 years running for both muon polarities, showing 
exclusion plots both in the {0[f^™\6i4)- and in the (^24, ^34)-planes. However, due to 
the impressive statistics achievable at the Neutrino Factory, it could well be possible 
that a positive signal is found (if sterile neutrinos with O(leV^) mass difference with 
respect to active ones do exist). For this reason, in Sec. 14.31 we have shown the region 
of the parameter space for which it is possible to distinguish the (3+l)-model from 
three-family oscillations. 

Eventually, we will present in this section a first analysis of the precision achievable in 
our setup in the simultaneous measurement of mixing angles and CP-violating phases. 
We first focus on ^24, ^34 and on the CP-violating phase ^3. Notice that each of the 
three possible CP-violating signals in a four-family model is related to a different 
Jarlskog invariant, proportional to a different combination of the mixing angles. The 
Jarlskog invariant that depends on sin ^3 is, in our parametrization, proportional to the 
combination sin 26*23524534 sin ^3, as it can be seen in eqs. fll6lll7p . A measurement of 
S3 is thus possible only if both 6*24 and 6*34 are simultaneously non-vanishing. We will 
thus show 99% CL contours in the (6*34, (53)-plane for particular input pairs (^34, ^3) for 
fixed non- vanishing values of ^24- 

The measurement of (^34, ^3) is achieved combining data from the z/^ disappearance 
channel and the discovery channel. This analysis, of course, does not pretend to 

be as exhaustive as those that have been presented in the framework of the three-family 
model. In particular, we will not address within a comprehensive approach the problem 
of degeneracies in four-family models. Notice that this problem, extremely severe in the 
three-family oscillation studies at the Neutrino Factory (see, for example, Refs. [69]ll08j 
and |122j ). is expected to be even more complicated in a four-neutrino model. In the 
particular case of the 53-dependent CP-violating signal, that can be extracted using the 
i^^i ^fi and Vr channels, we do expect to observe at least degeneracies due to the 
(6*34, 53)-correlation (the so-called "intrinsic degeneracies", [69]); those dependent on the 



wrong reconstruction of the sign of the atmospheric mass difference I ^m^i (known as 
"sign degeneracies", [70]); and, eventually, those dependent on a wrong reconstruction 
of the "atmospheric" mixing angle ^23 octant (known as "octant degeneracies", ^107]). 

The contours in the (^34, 53)-plane have been obtained as follows: we have first computed 
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At long baselines we are not sensitive to the sign of the SBL mass difference Am|^. 
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the expected number of events for z/^ — > and oscillations in the four- 

family model for particular choices of the relevant parameters, ^24 = ^24, 6*34 = ^34 and 
(53 = ^3. We have then computed the expected number of events in the (^34, 53)-plane 
for the same oscillation channels in the four-family model, varying ^34 G [0, 35°] and 
^3 G [0, 360°]. The Ax^ is then computed as follows: 




= L h^.(^24, ^34, h) - iV,(^24, ^34, ^3) (20) 



where the minimum of the is, trivially, obtained for ^34 = 6'34; b-^ = b-^. As before, 
J runs over the different signals: the disappearance and the —>■ Vr discovery 
channels data, divided into 10 energy bins, for the two baselines and the two possible 
stored muons polarities. The variance is defined by eq. ffT^ . with fj = 5% for the 
disappearance channel and 10% for the —>■ Vr discovery channel. No correlated 
systematic errors have been considered in the plots of this section. The region in the 
(6*34, ^3)-plane compatible with the input values (^34,^3) at the 2 d.o.f.'s 99% CL is 
eventually defined by drawing the contour line corresponding to Ax^ = 9.21. Notice 
that we have also studied the simultaneous measurement of ^^24, ^'34 and ^3 using the 
combination of the z/^ and Vj. channels, finding that sensitivity to ^3 is 

lost for values of the product S24S34 smaller than (s24S34)min ~ 0.01. The results that 
we show have been obtained for choices of the input parameters (^24, ^34) such that 
S24S34 > (S24S34) mill- 

We show in Fig.[12]the 2 d.o.f.'s 99 %CL contours for the simultaneous measurement of 
^34 and (^3 using the combined data from the disappearance and the discovery channels 
for two representative values of 6*24: ^24 = 3° (left panels) and ^24 = 5° (right panels). 
Blue dashed lines stand for the L = 3000 km baseline; red dashed lines stand for the 
L = 7500 km baseline; black dashed lines stand for the combination of both baselines. 



In the numerical analysis, the following parameters in the four-family model have been 
kept fixed to their central values: 6112 = 34°, 6113 = 0; Am'^^y = 7.9 x 10"^ eV^; Am|^ = 
Aml^^ = 2.4 X lO-^eV^; 5i = 82 = 0. Eventually, Am^^^ = 1 eV^ and Ou = 0°. For 
simplicity, we have fixed ^23 = 45°. We do not expect any "octant degeneracies", thus. 
The input values that we have studied to illustrate the discovery potential of our setup 
are: ^34 = 20°, 30°; 63 = 90° (upper panels) and ^3 = 200° (lower panels). Matter effects 
have been included considering, as always, a constant matter density p = 3.4 g/cm^ 
for the shortest baseline and p = 4.3 g/cm^ for the longest one, computed averaging 
over the density profile in the PREM [86] along the neutrino path. 

First of all, we can see that the combination of the two channels at the shortest baseline 
(blue lines) is not enough to solve the sign degeneracies (labeled with "SD" in the plot), 
that can be observed for all of the choices of the three input parameters (^24,^34,^3)- 
The sign clones are located at the point (6'f^,5f°), where 9f2 ~ ^34 and satisfies 
the relation sin sin A31L = — sin5|° sin A31L, with ~ —90° for ^3 = 90° and 
~ 20° for 63 = 200°. The intrinsic de generacy is also found for one specific choice 
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Fig. 12. 99% CL contours for the simultaneous measurement of 634 and 63 using the combined 
data from the disappearance and the — > Vr discovery channels. Two different values 
of 024: have been considered: ^24 = 3° (left panels); O24 = 5° (right panels). The input pairs 
(034, ^3 ), marked by a star in the plots, are: 634 = 20°, 30°; ^3 = 90° (upper panels) and 200° 
(lower panels). In the plots, "ID" stands for "Intrinsic Degeneracy"; "SD" stands for "Sign 
Degeneracy" . Blue dashed lines represent the L = 3000 km baseline data; red dashed lines the 
L = 7500 km baseline data; black dashed lines stand for the combination of both baselines. 

of the input parameter (^24 = 3°, ^34 = 20°, ^3 = 200°). On the other hand, no intrinsic 
or sign degeneracy are found at the longest basehne (red hues). When combining the 
two basehnes we see that the degeneracies are solved and that a very good precision 
on the simultaneous measurement of 634 and 63 is achieved for all the choices of the 
input parameters that we have considered. In particular, the error in ^3 at the 99% CL 
is of the order of a few tens of degrees. At the same time, the mixing angle 634 can be 
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measured for these particular inputs with a precision of a few degrees. 

We summarize our results for the simultaneous measurement of ^34 and ^3 in Fig. [T3l 
where the 99% CL " 53- discovery potential" in the (6*34, ^s) -plane for different values of 



^24 is shown J ^^1 We define the " (Js-discovery potential" as the region in the (sin^ 26*34, 63)- 
plane for which a given (non-zero) value of the CP-violating phase ^3 can be distin- 
guished at the 99% CL (for 2 d.o.f.'s) from the CP-conserving case, i.e., ^3 = 0, vr. Note 
that we have also taken into account the effects of the sign degeneracy in this analysis. 

In the left panel, only data from the u^j, —>■ disappearance channel are shown. In 
the right panel, we have combined data from the disappearance channel with those 
from the z/^ —>■ Vj. appearance channel. Upper panels refer to ^24 = 3°; lower panels to 
^24 = 5°. Blue dashed lines stand for the L = 3000 km baseline; red dashed lines stand 
for the L = 7500 km baseline; eventually, black dashed lines stand for the combination 
of the two baselines. 

We can see from Fig. [T3](left) that, using z/^ disappearance channel only, we are able 
to measure a non- vanishing ^3 for values of 6*34 above sin^ 2^34 > 0.4(^34 > 18°). The 
CP- cover ag is ~ 50%, with a very smooth dependence on ^34, being a bit larger 
for larger ^24- We can also see that the detector at L = 3000 km have no 53-sensitivity 
whatsoever. 

The situation is completely different when the z/^ — >■ discovery channel data are 
added to the z/^ disappearance ones. Fig. [131 right). First of all, we see that the L = 3000 
km detector is no longer useless to measure 53: spikes of (Js-sensitivity for particular 
values of ^3 can be observed, in some cases outperforming the far detector results. 
However, it is in the combination of the two baselines where we can see that a dramatic 
improvement in the 53-discovery potential is achievable. When the z/^ — > Uj. data are 
included, a non-vanishing ^3 can be measured for values of ^34 as small as sin^ 2^34 = 
0.06(^34 = 7°) for ^24 = 5° and sin^ 2^34 = 0.10(^34 = 9°) for ^24 = 3°. For sin^ 2^34 > 
0.4(^34 > 20°), roughly 80% (60%) of CP-coverage is achieved for ^24 = 5° (3°). The 
striking improvement in the 53-discovery potential is a consequence of the synergy of 
the two channels and of the two baselines, whose combination is able to solve most of 
the correlations that otherwise strongly limits the potential of the z/^ disappearance 
channel. 

For completeness, we also present in Fig. [H] results for the sensitivity of the golden and 
silver channels to the phase ^2, that reduces to the three-family CP- violating phase S 
in the limit 6*^4 ^ 0. A comment is in order: as it can be seen from eq. 0131) . using 
the parametrization in eq. ([2]), the golden channel oscillation probability in vacuum 
depends on the combination (^2 — ^3) up to the eighth-order in e. This means that, 
in the (3-|-l)-model, a CP-conserving result in the golden channel may be found for 



Notice that, as we stressed at the beginning of this section, sensitivity to ^3 is lost when 
the product S24S34 is smaller than (s24S34)min ~ 0.01. 

The CP-coverage is the fraction of the 53-parameter space for which we are able to exclude 
(^3 = 0, vr at the 99% CL for a given value of ^34. 
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Fig. 13. The 99 % CL "iJa -discovery potential" in the (634, dsj-plane. Left: only —>■ 
disappearance channel data; Right: combination of i^^ — disappearance and 
appearance channels data. Upper panels have been obtained for 624 = 3°; lower panels for 
024: = 5° . Blue dashed lines stand for L = 3000 km baseline data; red dashed lines stand for 
L = 7500 km baseline data; black dashed lines stand for the combination of the two baselines. 



non-vanishing values of 62 and ^3 if 61 = 0; (^2 — ^3) = 0, vr. This degeneracy can be 
broken only by adding new information, such as that obtained using the silver channel, 
see eq. f lT^ . disappearance or z/^ — > u^. appearance data. Golden channel data at the 
L = 3000 km baseline may not be able, thus, to detect a non-vanishing CP-violating 
signal even when both S2, S3 are different from 0, vr. 

The contours in the {9^^'^\ 62)-p\a,ne have been obtained as follows: we have first 
computed the expected number of events for z/g —>■ and z/g — > Vt oscillations in the 
four-family model for particular choices of the relevant parameters, Of^^™'' = ^(3^'^™'' and 
62 = 62- We have then computed the expected number of events in the {0[^^™\ 52)-plane 



for the same oscillation channels in the four-family model, varying ^if") G [0, lO'^ 
S2 G [0,360°]. The Ax^ is then computed as follows: 



and 




(4fam) 



S2)-N,{9 



(4fam) 
13 5 



(21) 



where the minimum of the is, trivially, obtained for 6 



(4fam) 
13 



e 



(4fam) _ 



13 



62 = 62- As 



before, j runs over the different signals: the z/g z/^ and the z^g data, divided into 
10 energy bins, for the two baselines and the two possible stored muons polarities. The 
variance aj is defined by eq. ffT^ . with fj = 2% for the golden channel and 10% for 
the silver channel. No correlated systematic errors have been considered in the plots 



of this section. The region in the 



3(4fam) 

corresponding to Ax^ 



T(4fam) 
'13 5 



52)-plane compatible with the input values 



{Oi^''"^' , 62) at the 2 d.o.f.'s 99% CL is eventually defined by drawing the contour line 

9.21. 
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Fig. 14. 99% CL contours for the simultaneous measurement of 613 and 62 using the combined 
data from the and the Vf, — > v.^ golden and silver channels. Two different choices of 

the active-sterile mixing angles have been considered: 614 = 624, = 5°; 634 = 20° (left panel); 
9iA = ^24 = ^34 = 10° (right panel). The input pairs (013,62), marked by a star in the 
plots, are: 613 = 2° , 5° ; ^2 = 90° , 250° . Golden and silver channels data are always summed. 
Four-family results are shown for the two baseline separately and summed: blue dashed lines 
stand for the L = 3000 km baseline data; red dashed lines stand for the L = 7500 km baseline 
data; black dashed lines stand for the combination of all data. Eventually, black solid lines 
stand for the three-family results for the combination of the baselines. 



Results have been obtained for 0^'^^^™^ 



'13 



90°, 250°. The other parameters 



2.4x10"^ eV^ and Am 
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1 



2°, 5° and 62 

are: 9^2 = 34°; ^23 = 45°; Am^^ = 7.9 x lO'^ eV^, An 
eV^. For simphcity, we show resuhs for 5i = ^3 = in these plots (remember that the 
measured phase should be interpreted as 6 = (62 — ^3) for 61 = 0). Eventually, the 
three active-sterile mixing angles are: 6*14 = 6*24 = 5°; 6*34 = 20° in Fig. [T^left): and 
^14 = ^24 = ^34 = 10° in Fig. [W right). In the plots, golden and silver channels data 
are always summed. Four-family results are shown for the two baselines separately and 
summed: blue dashed lines stand for the L = 3000 km baseline data; red dashed lines 
stand for the L = 7500 km baseline data; black dashed lines stand for the combination 
of all data. For comparison, black solid lines stand for the three-family results for the 
combination of the baselines. 

We can see from the plots in Fig. [14] (left) that, when the active-sterile mixing angles 
^14 and ^24 are "small", the four-family results are extremely similar to those obtained 



by a fit in the three-family modelj^ The 99% CL contours in the four-family model are 



slightly larger than the three-family ones. The shape of contours in the {Oi^^™\ (52)-plane 
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Results do not depend significantly on ^34: as it can be seen in eqs. ([T3]l and (fH]) . this 
angle only enters in the silver channel oscillation probability, statistically less relevant than 
the golden channel data. 
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is identical for the four- and three-family. This means that the correlation between the 
two parameters is not modified by contributions proportional to ^14, 624, with respect 
to three-family expressions. 

In Fig. [m (right) we see that, when 6^14 and 6*24 assume values near their upper bound, 
the results for the four-family contours can significantly differ from the three-family 
ones. In particular, the four-family contours in the (^[3^'^™'', 52)-plane are orthogonal 
to the three-family ones. This is easily understood by looking at the approximated 
expressions for the oscillation probabilities of the golden and silver channels expanded 
to order in Sec. 14. 1[ As we can see, the two terms with a (52-dependence in the 
second and third lines of eq. (fT3l) are proportional to S13S14S24 sin(52 + A31L/2) and 
to si3(A2i-L) cos(52 -|- A31L/2), respectively. When the first term becomes as large as 
the solar-suppressed one, the two terms give a destructive interference. Eventually, 
if the first term becomes larger than the second, the (^13^'^™'', 52)-correlation changes, 
becoming orthogonal to the three-family-like one, as it can be seen in the figure. Notice 
that, when 6^14 and 6^24 assume values near their upper bound, there is a little sensitivity 
to the CP-phases at the longer baseline. Since this baseline corresponds to the magic 
baseline, we can conclude that this is an effect characteristic to the four neutrino 
scheme. 



5 Conclusions 

We have studied the potential of a Neutrino Factory to search for signature of the 
(3+l)-scheme, where the largest mass squared difference Arngg^ can be any value larger 
than 0.1 eV^, as long as oscillations driven by Am^g^^ are averaged for the range of the 
energy and the baseline at a Neutrino Factory. From the analytic expressions of the 
oscillation probabilities, we have seen that the disappearance channel — > and 
the ^^discovery channeT —>■ Vr are extremely powerful in constraining the parameter 
space of the (3-|-l)-scheme. 

We have performed a numerical analysis of the sensitivity to sterile neutrinos in two 
setups: (a) A Neutrino Factory with the muon energy = 50 GeV, 2 x 10^° useful 
muon decays per year per baseline, with two detectors located at L = 3000 km and 
L = 7500 km from the source, respectively; (b) The setup suggested in the final ISS 
Physics Report, i.e., a Neutrino Factory with = 20 GeV, 5 x 10^° useful muon decays 
per year per baseline, with two detectors located at L = 4000 km and L = 7500 km 
from the source, respectively. The two detectors are of the Hybrid-MIND type (50 kton 
MIND -|- 4 kton MECC). In both cases we have assumed data taking for 4 years for 
both muon polarities. We have carefully taken into account the relevant backgrounds, 
efficiencies and systematic errors and our analysis is, thus, much more detailed than in 
previous studies in the literature. 

We have looked at the sensitivity to ^13, 614, analyzing the golden and silver channels. 
While Neutrino Factories do not give any useful constraint on 614 when marginalizing 
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over 024 and 634, as in the case of the three- flavor framework, they give quite a strong 
constraint on sin^ 2^13, down to 7 x 10"^ (2 x 10"^) for the 50 GeV (20 GeV) case 
with a shght dependence on ^14. In the present work we did not assume any near 
detector, since the issue has not yet been studied in detail (see |123j ). If we put a near 
detector, then we should be able to give a constraint on 6^14 because the near detector 
will give us information on how much neutrino oscillations occur due to the largest 
mass squared difference Am^gj^. Further studies on the near detector will be necessary 
to obtain quantitative results on 614^. 

Combining the disappearance and discovery channels, we have found that both the 50 
GeV (20 GeV) Neutrino Factories can constrain 6*34 and 6^24 down to 12° (14°) and 7.5° 
(8°) , respectively, marginalizing over all other four-family parameters. If we closely 
look at the details, then we found some difference between the performances of the 50 
GeV and 20 GeV Neutrino Factories, such as the sensitivity plot in the (6^24, 6'34)-plane 
or the regions in the (^24,^34)- and {0[^'^^\ ^i4)-planes where the three-flavor mixing 
hypothesis is excluded, where the sensitivity of the 50 GeV Neutrino Factory is always 
a bit better because the — >■ channel, which becomes more important at higher 
energies, makes this difference. The present results should be compared to the previous 
study [60] on the CNGS experiment, which cannot improve the present bound on ^24 
and ^34, if running at the nominal luminosity. The sensitivity to 6'24 and 6'34 of the 
Neutrino Factory setups proposed in this paper is better than the potential sensitivity 
of the CNGS experiment even if the CNGS beam intensity would be increased by a 
factor 10. 

We have also discussed dependence of the sensitivity of the four channels on the sys- 
tematic errors. The sensitivities of the golden and silver channels depend mainly on 
the correlated systematic error (Tq, on the overall normalization. As for the disappear- 
ance and discovery channels, the sensitivity of the former depends on the bin-to-bin 
uncorrelated error fj and cTq,, whereas that of the latter depends mainly on fj. To take 
advantage of the discovery channel, it is necessary to reduce fj below 10%. To this 
purpose, the expertise on EGG technology that will be gained after the first years of 
OPERA data taking will be of great importance. 

On the other hand, in the case where we have a positive signal of the sterile neutrino 
mixing, we have found that we can measure the new CP phase ^3 by combining the 
disappearance and discovery channels at the two baselines, where information from 
the L = 7500 km baseline plays a dominant role in resolving parameter degeneracy. 
After combination of the two channels and the two baselines, the CP-violating phase 
^3 can be measured at 99% CL for values of sin^ 2^34 > 0.06 (notice that, using only 
the disappearance channel at the two baseline, ^3 can be measured at 99% CL for 
sin^ 2^34 > 0.4, only). It should be emphasized that the measurement of CP- violation 
in the z/^ z/^ channel is a clear new signal of CP-violation associated to the sterile 
neutrino scheme. We have also found that the measurement of the three-family-like 
CP- violating phase 62 is modified by the presence of non-vanishing active-sterile mixing 
angles. 
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While our works contain some results on parameter degeneracy in the (3+l)-scheme, 
the problem on the general structure of parameter degeneracy in the four neutrino 
schemes is beyond the scope of the present paper and should be pursued in future. 

Finally, we would like to stress that, while the discovery channel at a Neutrino Factory 
is not very useful for the measurements of the three-flavor oscillation parameters, it is 
a very important channel to search for new physics beyond the standard scenario. The 
discovery channel at a Neutrino Factory deserves, thus, further studies. 



Appendix 



A The mixing matrix elements Uaj 

The mixing matrix elements in the parametrization ([2]) are given by the following: 



U, 



til 



u, 



/x2 



u, 



ti3 



u, 



m4 



-C23C24Sl2e 

C12C23C24 - Sue' 

Cl3C24S23e-*'^3 - 
C14S24 



Uel = C12C13C14 
— Ci3Ci4Si2e 

Uei = ■S14 



i(5i 



(A.l) 



C12 



C24Sl3S23e*^''2 ■'^^ + C13S14S24 

C24Sl3S23e*('''"'^3^ + C13S14S24 
162 



(A.2) 



■Sl3'Sl4"S24e 
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Us, 



si2e 

—C12 IC13C24S14S34 + Si3e 



C23C34 ~ ■523S24S34e 



-C12 



C34S23e*'^^ + C23S24"S34 



C13C24S14S34 + sise 



-C24Sl3Sl4S34e 
C14C24S34 



-iS2 



+ Ci3 



C23C34 - S23-S24S34e' 



C23C34 - 523^245346 



Sue 

~Ci2 IC13C24C34S14 — Si3e 



C23C34"S24 ~ 5235346 

1(52 



C34523524e + C23S3A } 



-C12 



C23C34524 - 5235346 



-Sue "-^^ {C13C24C34514 - 5i3e*'^2 



— C24C34Si3Si4e 
C14C24C34 



■iS2 



Cl3 



C34523524e + C23S3A 
C34523524e + C23534 



where c 



cos 9ij and Sij = sin 6ij. 



(A.3) 



(A.4) 



B Oscillation probabilities by the KTY formalism 



To derive the expressions for the oscillation probabilities in matter, we use the KTY 
formalism which has been introduced in Ref. The evolution equation of flavor 



eigenstateili 



is: 



.d 



where 



^ = diag(0, ^) ^ diag(0, A21, A3,, A4i^ 

^ = V2G'Fdiag(ne, 0, 0, 7x^/2) = diag(Ae, 0, 0, A^), 



:b.ii 



Aij = Am'fj/2E, and Ue and n„ are respectively the electron and neutron densities. 



In eq. ( IB. II) . we have subtracted from Ti. the term Eil = \J E"^ + mf 1, which con- 
tributes only to the phase of the oscillation amplitude and therefore does not affect the 



Another proof of the KTY formalism was given in Ref. |124|125j and it was extended to 
four neutrino schemes in Ref. [126J. 

Greek (Latin) indices label the flavor (mass) basis: a = e, //, r, s (i = 1, 2, 3, 4). 
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probability. In the KTY formalism, the oscillation probabilities in matter assume the 
following form: 



Pap = <^a/9 - 4 E Re{Xfxf*) sin^ (^^] + 2 ^ ImiXfxf*) sin( A^-L) , 



Kj 



\ 2 



(B.2) 



where AEji = Ej - E, and Xf = Uafipj (j = 1,2,3,4). Ei and Uai are the Ti 
eigenvalues and the effective mixing matrix in matter, respectively, defined through 

n = Udiaig(Ej)Ul 



The Xf^ matrices can be expressed as follows: 



where V is the Vandermonde matrix: 



i-i 



a(3 



(B.3) 



V 



/i 1 1 1 ^ 

El E2 E^ £■4 

Tp2 rp2 Tr2 t72 
^1 ^2 ^3 

\ElElElEl) 



whose determinant is Y{i<j ^Eji- The inverse of V can then be easily obtained as long 
as we know the eigenvalues Ej of the effective Hamiltonian in matter Ti, expressed in 
terms of Ag, An and the vacuum parameters: 



/ 



1 



Aii/21 Ail/41 

A£'2l A£'32 A£'42 
AE'si Ai^32Ail^43 
V A^4lA^42A^43 



(£^2-£'3-E'4, —{E2E3 + E^E^ + E4E2), E2 + E^ + E4, —1) 
[E3E4E1, —{E3E4 + E4E1 + E1E3), E3 + E4 + El, —1) 
(E4E1E2, —{E4E1 + E1E2 + E2E4), E4 + El + E2, —1) 

{E1E2E3, —{E1E2 + -E'2-E'3 + E^Ei), Ei + E2 + -Es, —1) 



(B.4) 



Within the KTY formalism, thus, we only need to compute the eigenvalues of 7i to 
derive the oscillation probabilities in matter. 

A possible drawback of this approach is that the physical understanding of the oscil- 
lation probabilities (i.e. the dependence on the vacuum mixing matrix parameters) is 
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encoded in the explicit expressions for the X coefficients. To make contact with the 
parameters to be measured in a manageable way, we thus need to introduce some ap- 
proximations in the computation of the eigenvalues Ei and of the corresponding matri- 
ces X°'^ . Now, considering the present constraints from |60j in the standard and sterile 
small parameters, we see that ^13,^14 and 624 cannot be much larger than 10° while 
the third active-sterile mixing angle, 634,, can be as large as 634 ~ 35°. Notice also that 
the present constraint on the 623 deviation from the maximal mixing, 6623 = ^23 — 7r/4, 
is of the same order as those on 6*13, 6*14 and 6*24. On the other hand, the solar and 
atmospheric mass differences, Arn^^i, Aml^^^, are much smaller than Am^^-^. In what 
follows, therefore, we expand all the quantities in power of a small parameter e, and 
keep terms of cubic order in e, where the small parameter is defined by 



e = 6'34 ~ \/6'i3 ~ y 6*14 ~ Y 6'24 ~ V 56*23 < 4 X 10 \ 

T]2 = AmlJAml,<10~\ 
r]3~AmlJAml, < 10'^ 

Notice that, to third order in e, in the expansion in the probabilities we have neglected 
all terms proportional to ?7e,n,2,3- Although this can be a rather rough approximation, 
as we have seen before, it is very useful in order to understand the different physics 
potential of the various oscillation channels. Thus we have the following probabilities 
to third order in e: 



Pee~l + 0(e^), 

Pe,^Per ~ Pes ~ O (e^) , 

-P/iAt = 1 - sin^ - 2 (AnL) S24 S34 cos ^3 sin A31L + O (e"^) , 

PfiT = (1 - -534) sin^ + {s24 S34 sin ^3 + 2 (v4„L) S24 S34 cos 63} sin A31L 



P/is = S34 sin^ - S24 S34 sin 63 sin A31L + O (e^^ 



In Sec. I4.1l we had to go beyond O(e^) in order to explain our numerical results using the 
golden and silver channels. To this purpose, in the text we have shown the approximated 
expressions for Pg/^ and Per to order in vacuum. To check unitarity of the four-family 
PMNS matrix at this order in e, it is useful to show here the Pes oscillation probability, 
also: 
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Pes — 2 6*14(1 - 6*14 - 6*24 - ^34) 

+ 2 {^^^^(-2 el + + - 25^23^14) + ^?3^24^^34 COS ^3} sill 



2 Asii^ 



2v^ ^13^14^24(1 + Se23 - 0I4) Sin ( ^2 - ^3 + ^|^) sin 



2x/2 ^13^14^34 - 5^23 - + 



A,iL 



A31L 



31J^ \ . ^A3iJ 

sm ■ 



- sin 26^12 6'i36'|4(A2iL) cos (^6i -62 + S3 - sin 

+ sin 2^12 ^i4^34(A2i-L) sin(5i + '^s) - sin 2^12 ^i4^24(A2i-L) sin 5i . 

As a final analj^ical contribution, we have calculated approximate probabilities asso- 
ciated to the channels under study, P^^ and P^,- (together with Pf^s), to fourth order 
in e but neglecting ^13 and ^14: 



-P^i^i — 1 2 ^24 



1 - 4{se2sr - + ol.^ (4<^^23 - ^ 

^31 V 



34 A 
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